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Absoraco 


Undesired gravitational affects scan as convection or sedimentation in 
a fluid can sometimes be avoided or decreased by tr.e use of a closed cnamoar 
uniformly rotated about a horizontal axis. 

In a previous study (Roberts, Kcmield and rcwlis, 1939) we determined 
the spiral crtits of a heavy or buoyant particle in a uniformly rotating 
fluid. The particles move in circles, and spiral in or cut under the combined 
effects of the centrifugal force and centrifugal buoyancy. We also formulates, 
and solved an optimization problem for the rotation rate of a cylindrical 
reactor rotated about its axis and containing distributed particles. 

This report is concerned with related studies in several areas. 

We have upgraded a computer program based on our analysis, correcting 
some minor errors, adding a sophisticated screen- and-pr inter graphics 
capability and other output options, and improving the automation. 

We have supported the design, performance, and analysis of a series of 
experiments with monodisperse polystyrene latex micrcspheres in water, to test 
the theory and its limitations. 

The theory was amply confirmed at high rotation rates. But at low 
rotation rates (lrpm or less) the assumption of uniform solid-body rotation of 
the fluid became invalid, and there were increasingly strong secondary motions 
driven by variations in the mean fluid density due to variations in the 
particle concentration. In these tests the increase in the mean fluid density 
due to the particles was of order 0.015%. 

To a first approximation, these flews are driven by the buoyancy in a 
thin crescent-shaped depleted layer on the descending side of the rotating 
reactor. This buoyancy distribution is balanced by viscosity near the walls, 
and by the Coriolis force in the interior. A full analysis is beyond the 
scope of this study. 

Secondary flews are likely to be stronger for buoyant particles, which 
spiral in towards the neutral point near the rotation axis under the influence 
of their centrifugal buoyancy. This is because the depleted layer is thicker, 
and extends all the way around the reactor. 
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Section 1 
TNTRC CUCTICN 


1.1 Background 

The sffecus of caroicie or ceil sedimentation m a fluid can ce 
minimized by rotating the fluid container about a horizontal axis. The 
container itself need net necessarily be axisymmetric. 

been used for experiments and for materials 
alone, and as a preliminary to micrcgrav ity 

If the fluid rotates uniformly with its container, then individual 
particles mcve in circular orbits, with the same rotation rata, about a canter 
displaced horizontally from the axis cf rotation. For particles at this 
center, their sedimentation velocity just cancels the upwaro or downward ficw 
of the surrounding fluid. For particles denser than the fluid, the orbit 
radius increases slcwly with time due to the net centrizugal force, so that 
the particles spiral outward from their center. Buoyant particles spiral 
inwards, due to the net buoyancy force. 

Relative to rotating axes fixed to the container and fluid, all the 
particles with a given sedimentation rata move in phase, in circular orbits 
with a fixed ra dius equal to the separation between the axis and the center 
described above. The velocity of every particle is its sedimentation rate, in 
the current direction of gravity. This circular motion is superposed on a 
centrifugal motion away from or towards the axis. Particles with larger 
sedimentation rates move in larger circles. 


The technique has 
processing. It is useful 
exceriments in space. 


1.2 Applications 


others: 


The technique has been used in the following applications, among 


Production of latex microspheres of uniform size; 

Study microgravity effects on cells and tissue; and 

Minimize gravitational flow disturbances in free-flow 
isoelectric focusing. 

Other potent ial applications include crystal growth (where convection is 
sometimes a problem) . 
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The alternatives of micrcgravity processing in space flight, in a 
plane flcwn on a parabolic flight path, or in a drop tower, are generally 
either too inconvenient and expensive, or too brief. 


1.2 Requirements 

Successful use of this technique in regard to particle processing (as 
for latex) imposes the following requirements: 

Keep particles from collecting at container boundary before 

processing; 

Keep particles from collecting at container boundary during 

processing; 


Selective withdrawal of particles after processing, to exclude 
these that have interacted (or interacted excessively) with the 
container boundaries; and 

Temcerature control, where calls, tissues, or processes require 
it. 

In addition, it is necessary to determine and increase any upper limits on the 
particle concentration that can be processed without seriously modifying the 
flew. 


1.4 Subjects for Analysis 

Experimental and theoretical study is needed in regard to the 
following questions, among others: 

What happens to particles that hit the container boundary? 

Are non-axisymmetric containers useful? For cylinders, what 
size and aspect ratio is appropriate? 

Can useful flows be established during initialization by 

changes in the rotation rate? These spin-up and spin-down 
effects are dominated by the Coriolis force. 

How is the mean flow modified by the aggregate net weight of 
the particle distribution? At what stage does this 

modification became a problem? An example of a related study 
is the work of Batchelor and Janse Van Rensburg on bidisperse 
sedimentation (1986) . 
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1.5 Methods of Study 

These problems should be studied using an appropriate liCmbinaticn of 
experiments , analytic methods, and computer models. A cicse interaction 
between these approaches is appropriate. 

Analytic methods include solving ordinary dif-sreni-ial equacicn 
sys tem s ( as for orbits, cf. Roberts eu al. 1989) , or partial differential 
equations (linearized calculations of rotating flows usmc a balance cf 
buoyancy, viscosity forces, and pressure gradients). 

Ccmcuter models are generally required when the prcolers are too h ar a 
for analytic methods. Cur Atmospheric General Circulation Experiment (AGv_i) 
computer code, developed under MSFC sponsorship, has great flexibility in the 
study of a wide range of rotating flews, with or without temperature 
variation, and has previously been applied to spin-up and to the isoelectric 
focusing cas e described above. In the liiruxed scope of this program, computer 
modeling was limited to orbit calculations. 

Frequent interaction with laboratory workers and involvement in the 
experiments should be a part of every such study. Cur practice has been to 
maintain a close relationship with the laboratory workers, and to support the 
design of the equipment, the choice cf the experimental procedures and 
techniques, and the interpretation of the observations and results. 
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Section 2 


RESULTS 


2.1 Experiments with Polvstr/rene Latex Micro spheres 

We have supported the design, performance, and analysis of a series of 
experiments on reactor performance. Thes experiments were dene by Dale 
Kcrnfald at MSEC. Cur ccncrihuticns included assisting in the improvement of 
the optical systems, and in setup and performance of the difierant tests. 

The experiments were done with suspensions in water of mcncdisperse 
polystyrene latex micrcspheres, with diameters from 30 to 100 microns. The 
particle concentrations were of order 0.1% by volume. The reactor was a 
cylinder, with interior diameter 2.125 inches and length 3 inches. The curved 
boundary and one plane end boundary were transparent. The non- transparent end 
of the bath was a piston, moving in and cut to compensate rcr changes in the 
fluid volume in the cylinder. 

The r ea ctor was set up with water, and allowed to reach uniform 
rotation. High-concentration latex suspension was then injected with a 
syringe through the moving piston. The syringe could be operated in and cut 
with sufficient force to flush the syringe and to create turrulence inside the 
chamber, with thorough mixing. Within a minute of in^ecoicn, observation 
shewed uniform particle concentration, with no apparent deviations from solid 
body rotation of the fluid. 

The particles were illuminated in a darkened laboratory by a vertical 
slice of laser light with thic3<ness about a millimeter, entering horizontally 
and normal to the chamber axis. Refraction problems were minimized by filling 
the surrounding water bath. The particles in the illuminated slice could be 
viewed directly through the sides of the cylinder, and through a rotoscope at 
the end of the cylinder (which eliminated the rotation and made the chamber 
appear stationary) . The light was strongest at small scaooering angles from 
the incident light. The most useful results were obtained using a video 
camera, with a special cylindrical lens to make the image of the illuminated 
slice appear circular even though it was viewed at a small scattering angle. 
This led to focusing problems; our only solution to date has been to use so 
much light that the camera aperture was very small. 

Until now, all our experiments have been done with latex particles of 
density 1.05, heavier than the fluid which was water. We are considering 
using the same latex particles in a fluid with a density greater than 1.05. 

The first group of experiments were designed to confirm the rate of 
centrifugal spiralling outward, and therefore used relatively large rotation 
rates, of order 20rpm. At these rates, no secondary flow. effects were 
observed. The exponential decrease in the particle concentration with tune 
was measured, and was in agreement with the centrifugal settling theory. 
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The second gro up of experiments 'used lower rotation rates, from about 
0. 6rpm to 2rpim. At these rates, the centrifugal spiralling outward is 
relatively minor, and the depleted crescent formation is in theory likely to 
be more significant. This was confirmed, with a visible crescent on the 
descending side at 1.5rpm. 

But at lower rotation rates, sarong secondary flew effects became 
accarant, with time and space variations in the particle concentration and 
motion. The change with decreasing rotation rate was quite abrupt, especially 
at cur largest concentration of 0.3%. At 1 . 5rpm nothing interesting happened 
except for the crescent, while at 0.7rpn me secondary ficws grew rapidlv,^ in 
about 10 minutes, and moved most of the particles to the walls, leaving thin 
plumes of high concentration ascending at the ascending wall, and descending 
across the chamber interior under their own weight. The patterns varied on a 
time scale of about 2 minutes. large regions of the visible slice became 
almost void of particles. The results at lrpm shewed significant secondary 
flows, but they were not nearly so chaotic as those at 0.7rpm. 

The flow fields were three-dimensional; this is the only explanation 
for seme abrupt changes in the local concentration which were not advected in 
from adjacent parts of the visible slice. And they were driven by variations 
in the particle concentration. 


2.2 Theory of the Secon dary Flows 

We have analyzed the video res'll ts extensively . Particles with 
specific density 1.05, at 0.3% concentration by volume, increase the density 
of the water by about 0.015%. Even at particle concentrations as lew as 
0.03%, the nonuniformities in the particle concentration clearly drive a 
secondary flow, superposed on the uniform rotation, for rotation rates of lrpm 
and smaller. The nonuniformities in the particle concentration can arise only 
from the wall, since particle concentration is conserved following a fluid 
element, while decreasing exponentially at a low rate due to centrifugal 
buoyancy. 

From observation, once particles collide with the wall (or with the 
aggregate of particles against the wall) , only a small proportion are ever 
resuspended. The secondary flews are driven primarily by the net buoyancy of 
the depleted crescent, and secondarily by nonuniformities in the particle 
concentration associated with resuspensicn. Theoretically, these flows are 
proportional to the square of the thickness of the depleted crescent, which is 
in turn inversely proportional to the rotation rate. If these secondary flows 
become comparable with the solid body rotation, the theory based on uniform 
rotation breaks down completely, and the flows and particle distributions are 
quite different (and plainly three-dimensional in the videos) . 

The experiments make it plain that the lower limit on the rotation 
rate, predicted in the earlier theory based on the approximation of solid-body 
rotation, is not realistic for a concentration by volume of even 0.01%. But 
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concentrations of up 
mcncmer-swollen latex 
factor of order 100, 
diameter. 


to 30% are 'used in typical recipes for polymer isuticn of 
micrcscneres . The rotation rate must be larger by a 
so that the depletion layer thickness is only 1% of the 


This statement assumes that the reactor has a lerctr. comparable 
its diameter. If the length is very small compared with the ciameter, 
the secondary flews will be limited by viscous interaction wi_t the er.es. 


* 


h 


then 


In a latex reaction, the initial moncmer-swcllen particles have a 
relative density of order 0.9, ar.d tend to centriruge inwarts, away -om — e 
wails. Thev increase in density as toe absorbed Icv-aensity monomer is 
polymerized *to high-density polymer, according to a non- linear conversion 
versus time curve. This buoyancy distribution will also estaonsh its own 
secondary flew, and the overall result is hard to predict. 


2.3 Computer Solution for Uniform Rotation 

We have completed the computer program upgrade. 

We reviewed two computer programs written by otters. The f — s^, 
FAST. FOR, was written by a summer student employee under the supervision of 
William Fowlis, during 1987. It implements the analysis cf cur report, and 
for given keyboard input parameters computes the various resulting parameters 
as described in the report. It also writes output files in a suitable form 
for the commercial EnerGraphics program, which reads these flies and produces 
screen and printer plots of the spiral particle tracks. 

The second program was written by John Clelana, at the Research 
Triangle Institute, NC. It is identified as KORN2.PCR, named for Dale 
Komfeld, who manages this effort at MSFC. This program is a similar 
implementation of the analysis of our report, reading input parameters from 
the keyboard, and computing the various resulting parameters. It a^so 
generates simple plots of particles distributed at the computed concentration, 
at s uccess ive times, using the commercial Grapher software. 

We chose to upgrade the second computer program and to add superior 
graphics. Basically, the two codes perform similar functions. . The spiral 
plots produced by FAST. FOR are not particularly helpful, since in practical 
cases the thickness of the lines fills the figure. Similarly, it takes two of 
the K0RN2.F0R particle plots to reveal a single concentration change factor. 
And the use of the EnerGraphics and Grapher packages, in two stages, was 
unnecessarily complex. 

We corrected two errors in the KDRN2.P0R program. First, there was a 
minor error in the computation of the viscosity of water at lew temperatures, 
possibly by an error in transmission of the code. This is documented 
in Appendix A. Secondly, the approximate solution of a cubic equation for the 
parameter as a function of was seriously in error for larger values. We 
implemented a high-accuracy iterative solution, as described in Appendix B. 
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We ungraded the data input for the program. It new reads its 
parameters from a file containing the values used last, arc the user need only 
enter changes from these values . Each program execution will . compute cases 
one after the other, unci! the series is terminated cy entering a negative 
tire. Also, we implemented a variable reacccr radius. 

We have added optional graphic p^cts of tne restate. These grapnics 
are plotted using the ccmmercial software paexage Grapmu, from Sciett: — ic 
Endeavors Corporation, of Kingston, TN. For each case, two plots can be 
produced. 

The first graphic is illustrated in Figure 1. It displays a cross 
section of the rotating reactor,, with the fluid in sclid-ocdy rotation about 
the center 0. There are two other concentric circles, centered cn the point A 
where a stationary particle is just supported against its weight (net of the 
weight of the displaced fluid) by the viscous drag of the fluid flew past it. 
The interpretation of this figure is different depending cn whether the 
particles are heavier or lighter than the fluid. 

For particles heavier than the fluid, the fluid flew is 
counterclockwise in the figure, and is assumed to be solid— body rotation as 
described earlier. Particles at the point A remain stationary; with the 
downward force of gravity balanced by the upward buoyancy and drag forces. All 
other particles move outward from A in counterclockwise spiral orbits, with 
the ra di us increasing by a fixed and very small fraction on each orbit. 

Particles which start in the crescent shape hit the lever semicircle 
of the wall during their first orbit. Ideally they can fall off the wall in 
the upper right of the figure, but observation of latex particle experiments 
suggests that few do, and the crescent is cleared of particles in ten or so 
rotations . 

Heavy particles starting cut between the two concentric circles spiral 
outwards during the fixed time period of interest until they cross the outer 
circle and hit the wall. As stated earlier, after one or more such impacts, 
they adhere to the wall, and do net fall away. 

Heavy particles starting out inside the inner circle also spiral 
outwards, until at the end of the time of interest they fill the outer of the 
concentric circles. The particle concentration remains uniform during this 
process. 


For particles lighter than the fluid, the fluid flow is clockwise in 
the figure, and is assumed to be solid-body rotation as described earlier. 
Particles at the point A remain stationary; with the upward buoyancy force due 
to the weight of the fluid displaced balanced by the dewnward gravity and drag 
forces. All other particles move inwards towards A in clockwise spiral orbits, 
with the radius decreasing by a fixed and very small fraction on each orbit. 


- 3 - 


Rotating Reactor Studies 


Heberts Associates 


Buoyant particles which start in the crescent shape hit the uppe- 
semicircle of the wall during their first emit. In theory they can 1^-ft or 
the wall in the lower righr cf the figure, but the observations of lata, 
particle experiments (with heavy particles, as described abeve; suggest t 
St of then will adhere to the 'wall, with the crescent cleared of particles 
in tan or fewer rotations. 

Liont particles starting cut with unifora concentration inside the 
outermost cf the two concentric circles spiral inwards during the fixed time 
period of interest. By the end of the tire interval, they are ail inside the 
inner circle, with their concentration still uniform. 


The second graphic produced for each case is illustrated in Figura _ ^ • 
It shews plots cf the following three diners ionless functions of the rotation 

rate _Qj(in rpm) : 

r^/r Q = expC-OMi^) 2 ] = ex?(-6) ■ 3rd 

m/b = expC-OUl.^) 2 ] = (l-<5) exp(-6) - JP . 

Here r. and r are the radii of the inner and outer concentric circles, and b 
is the 1 reactor radius. The other notation is as in the report in Appendix D. 
We have chosen to maximize F, which is the ratio of the area of the snia.. 
circle to the cross section area of the reactor. Hius the optimum rotation 
rate Si. corresponds to the peak value of the third function. It should be 
noted that the shape of the first two functions is always the same, but tie 
form of their product depends on the relative values of dL^ and db^. 


2.4 Report on Uniform-Rotation Theory 

We have made changes to our earlier report "Particle Orbits in a 
Rotating Liquid" on rotating reactor theory, as requested by the referees, and 
resubmitted it to the Journal of Fluid Mechanics. The report is attached as 
Appendix D. This paper assumes solid-body fluid rotation, and neglects all 
secondary flows. It is the basis for the computer program described above. 
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Hours 

10.0 

Microns 

200.0 

RPM 

2.0 

Reactor Radius 

3.0 

Density Excess 

-0. 100 

Viscosity (p) 

0.0 LOO 

RESULT 

S 

| 

Area Fraction 

0.2 l 16 

Spiral Offset 

-1.039a! 

Radius Growth 

0.7040 ! 



Figure 1. Reactor Cross Section with the Inner and Outer Circles 

Particles starting in the crescent hit the wall immediately. 

Heavy pa rticles starting between the circles spiral out to hit the wall. 
Heav y particles starting in the inner circle spiral out to fill the outer 
Light particles starting in the outer circle spiral in to fill the inner. 
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DATA 


! Hours 10.0 

| Microns 200.0 

! J j Reactor Radius 3.0 

Density Excess -0.100 

i.T Viscositv 0.0100 


RESULTS 


Minimum 

RPM 

0.073 

Maximum 

RPM 

0.354 

| 1 Optimum 

RPM 

0.194 

! Cotimum 

F raction 

0.463 
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solution "can be creamed using the methods of Ekman layer analysis, ^mciucing 
end effects ac the cylinder ends. These solutions can ce confined and 
refined numerically, and with appropriate measurements. Such a program would 
provide a basis for optimizing the rotating reactor system and extending its 
range of application. 
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APPENDIX A 

\jj c CDS TTY FCPffiJLATTCN FOR WATER 


v^SCCSL~ 

increases frcm fr=e. 
processes 


cr wars 


nc 

Ox. cax 

prcducaicn cf ircncciispsi 
70 decrees C. 


to tci-ing. 
inoarsso ars 


:y a faces r or o as or. 
Ls criar.ce is 
re as. lew tairperatra 
=r*3s recuares rear era 


see tame 



reacs 


the oarxeracara froio ore 


syfccard , 


arc 


Cur ccixucar prccran - ~ — — _ , , 

commutes the visocsif/ Gsirg t-c empirical fcrmulae frem Seast ; 1 1989) .based 
on data from the National Eureau cf Standards. pe^ j-crtu-oS a P^ i i w 
temperature above and below 20 decrees C, and were inciuced i-~ '- he pregram 
KORN2 . FOR which we have upgraded. 

In validating the formulas, we 'discovered an error in their code 
implementation for lew temperatures . The following program _ and liszu^ 
computes the bad formula and the two good ones over tie ',,.uie range of 
temperatures, and lists also the tabulated values from Weas. (1S8S, . 

The following two points should be noted. First, tne. ccirputad 
viscosity was too large by only about 1% at zero degrees an msignir leant 
error. Secondly, the formulas agree with the taculat— _ values ^ for 

temperatures below and above 20 degrees, except that the tabulatec values were 
truncated to four significant figures above 20 degrees, instead of being 
rounded. 
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N33YI5C.FCR 


Fr ; c 3 v, Jan^arv 11, 1991 3 


0 pm 


program nbsvisc 

NATIONAL 3URE.AU OF STDS. FORMULAE FOR H20 VISCOSITY v -OISE, 
CR 1 GM/SEC CM); CRC HANDBOOK OF CHEM. AND =- v S., 53'- ED., 
pp . F - *9 TO F-51 . 


wri te(0,2) r , . 

formate COMPARISON OF NBS FORMULAS FOR VISCOSITY OF ,«-cR // 
1T3, 'TEMP 1 , T13, 'VI SCO ' , T25, 'VI SCI 1 , T37, ' V I SC2 ' , TA6, 'VI Si. , - VISC2'/) 
do 1 temp=0. , 100 . r 2. 


c 


c 

c 

c 

c 


c 


c 

1 


if ( temp. le. 20. ) then 

v=(1301 ./ (998. 333+8. 1 855* ( temp- 20. 30585* 

1 (( temp-20 . )*2) ) ) - 3.30233 

vi sc0=(10.0)**v 

PROGRAM KORN2. FOR AS RECEIVED DID NOT SQUARE THE SECOND TERM 
IN THE POWER SERIES 

v=(1301 ./(998. 333+8. 1855*< temp-20. )+0. 30585* 

1 ( ( temp- 20 . )**2) ) ) * 3.30233 

vi scl =( 1 0 . 0)**v 

else 

v=(1 .3272*( 20. -temp) -0.001 053* ( (temp- 20. )**2))/ 
1 (temp+105.) 

visc2=. 01002 * ( 10.0)**v 

endi f 

wri te(0, ' (f8.0,4f 12.3} 1 ) temp, vi scO , vi scl , vi sc2, vi sc 1 - vi sc2 


end 


CRiG:.NnL IS 

of poo» ouaittv 
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:o j c;riscn of mbs formulas for viscosity of 


:evp 

VISC0 

VI SCI 

VISC2 

j . 

.01306736 

.01736397 

.0 1776873 

2 . 

.0163596 6 

.01671371 

. 0 * 663208 


.01577390 

.01566359 

.01560313 

o . 

.01479456 

.0147204 7 

.01-66908 

3 . 

.01390846 

.01385306 

.01331371 

1C . 

.01310437 

.01307162 

.0*304250 

* 2 . 

.01237263 

.01235273 

.01233220 


.01170497 

.01169405 

.0**68064 

'6 . 

.01109422 

.01108921 

.01*03161 


.01053416 

.01053263 

.0105 2965 

22 . 

.01001941 

.01001941 

.01002000 

22 . 

.00954526 

.00954526 

. 00954847 

2-. 

.00910757 

.00910637 

.00911137 

26. 

.00870273 

.00869940 

.00870544 

28. 

.00832754 

.00832136 

.00832778 

30. 

.00797920 

.00796963 

.00797584 

32. 

.00765520 

.00764184 

.00764733 

34. 

.00735334 

.00733590 

.00734020 

36 . 

.00707162 

.00704992 

.00705264 

3c . 

.00630831 

.00678223 

.00678300 

mC . 

.00656182 

.00653130 

.00652981 

42 . 

.00633075 

.00629578 

.00629176 

4** . 

.00611383 

.00607444 

.00606764 

46. 

.00590991 

.00586617 

.00585638 

48. 

.00571797 

.00566995 

.00565699 

50. 

.00553708 

.00548488 

.00546860 

52. 

.00536640 

.00531013 

.00529039 

54. 

.00520516 

.00514495 

.00512163 

56. 

.00505267 

.00498865 

.00496165 

53. 

.00490830 

.00484060 

.00430984 

60. 

.00477148 

.00470023 

.00466564 

62. 

.00464167 

.00456702 

.00452854 

64. 

.00451840 

.00444048 

.00439807 

66. 

.00440123 

.00432016 

.00427379 

68. 

.00428975 

.00420568 

.00415531 

70. 

.00418359 

.00409665 

.0040422 7 

72. 

.00408241 

.00399272 

.00393432 

74. 

.00398590 

.00389359 

.00383115 

76. 

.00389377 

.00379895 

.00373248 

78. 

.00380574 

.00370854 

.00363805 

80. 

.00372158 

.00362211 

.00354759 

82. 

.00364105 

.00353941 

.00346089 

84. 

.00356394 

.00346024 

.00337773 

86. 

.00349006 

.00338439 

.00329791 

88. 

.00341922 

.00331168 

.00322125 

90. 

.00335124 

.00324193 

.00314758 

92. 

.00328598 

.00317498 

.00307674 

94. 

.00322329 

.00311069 

.0030085 7 

96. 

.00316302 

.00304890 

.00294294 

98. 

.00310506 

.00298949 

.00287972 

100. 

.00304927 

.00293233 

.00281878 


■.A , z x 

Vise - - 7 1 SC2 ~ABl3 IN EDI’ICN 'C , PAGE F- 40 

.00010024 var 
.;:gc3165 
. JG0C6541 
.00005133 
. 00003935 
. GO 0029 12 ^307 

.00002053 

.00001340 

. GCGG076Q 
. GGGQ0298 
- .GCG00059 1GG2 
- . SCO 0032 1 
- . 0CQGO5OO 
- .30000604 
- . G0000642 
-.30000622 7975 
-.30000549 
-.30000431 
-.30000272 
-.00000077 
.00000149 6529 
.00000402 
.00000680 
.00000979 
.00001296 
.00001629 5468 
.00001974 
.00002332 
.00002700 
.00003076 
.00003459 4665 
.00003843 
.00004241 
.00004638 
.00005037 
.00005438 4042 
.00005841 
.00006244 
.00006647 
.00007050 
.00007451 
.00007852 
.00008251 
.00008648 
.00009043 
.00009435 
.00009825 
.00010212 
.00010596 
.00010977 
.00011355 


VISC0 USES THE BAD T<20 FORMULA IN MY COPY OF KORN2. FOR (OMITTED * IN r-20)*’2) 
THIS MAY HAVE BEEN A TRANSMISSION ERROR. 

VISC1 IS CORRECTED, AND AGREES WITH THE 4 SIGNIFCANT FIGURES IN THE REFERENCE. 

VISC2 (FOR T>20) AGREES WITH THE REFERENCE, EXCEPT THAT THE REFERENCE TABLE 
VALUES ARE TRUNCATED TO 4 FIGURES INSTEAD OF ROUNDED! 
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APPENDIX 3 

SOHJTICN OF THE CUBIC ECGVTTCN FOR DELTA 


The two programs FAST. FCR and KCi 
urcr ads did net attempt to solve a cubic £ 
widen arises in cur analysis in the repero 
of one rotation rata. Physically ar.d mama 
cr.e as O increases from zero to iniinity. 
aepreximate formulations for ci, nominally v 
have replaced these with an erricient i 
crccram validates cur solution aigcrio 
approximations . 


RN2.FCR made available to us for 
muaticn for a as a function cf 
in Appendiuc A of the optimization 
raticallv, a increases from zero to 
The second program introduced^ three 
alia in different ranges of <t>. We 
terative solution. The following 
hn, and compares it with the 


The program takes a representative set of exact <T values, and for eac., 
cr.e it computes the corresponding Without using the known c, it computes 
the three empirical approximations and the first five iterations. The result, 
list the exact value, the empirical approximations, their errors, the <p Vai-e, 
the fifth iteration value, and the errors of the first five iteratm/e values. 


The following points should be noted. First, the empirical formulas 
are reasonable for <fc values up to about 1.3, with errors in d^of less than 
10%. Secondly, the first formula is better than the second for c values up to 
0.15, rather than the 0.05 in the supplied code. Thirdly, _ the empirical 
formula is a disaster for larger <£> values, predicting an unphysical 6 value of 
more than unity, and leading to program failure when a fractional power or 1 <3 
is e^/aluated. 
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CCWAR:S - N 0F CLHWNO "CURVE FITTING" SOLUTIONS FOR DELTA WITH EXACT VALUE AND SIMPLE NEWTON 


TERATICN SOLUTION 


DEL 

DELC1 

. OOGCO 

.00000 

.31 OCC 

.01003 

. 22C0C 

.02014 

■T 7 f, 

.03031 

^ n •"* 

.04055 

: c 

.05086 

.062 CC 

.06125 

.37000 

.07171 

. uS vWLwO 

.08225 

.09000 

.09237 

. 1 0 0 C 0 

.10357 

. *5000 

.15335 

.20000 

.21544 

. 25CC0 

.27516 

,30000 

.33787 

,35000 

.40405 

. 4C0QQ 

.47425 

.45000 

.54924 

50000 

.62996 

55000 

.71773 

60000 

.81433 

65000 

.92234 


.70000 

.75000 

.30000 

.35000 

.90000 

.91000 

.92000 

.93000 

.94000 

.95000 

.96000 

.97000 

.98000 

.99000 

.99100 

.99200 

.99300 

.99400 

.99500 

.99600 

.99700 

.99800 

.99900 


1.04566 
1.19055 
1 .36798 
1.59976 
.93899 
.03061 
.13513 
25658 
40110 
.57870 
.30706 
.12175 
.61035 
.59517 
. 76423 
.96000 
,19099 
,47019 
31880 
27441 
91281 
92113 
99000 


.99910 10.34812 
.99920 10.76356 
.99930 11.25460 
.99940 11.34920 10. 
.99950 12.59291 11. 
.99960 13.56666 11. 
.99970 14.93353 13. 
.99980 17.09634 14. 
.99990 21.54219 18. 


DELC2 

.00000 
.00923 
.01343 
. 02767 
.03694 
. 04626 
.05564 
.06507 
. 07457 
.03412 
.09374 
.14290 
. 19400 
.24735 
.30329 
.36224 
.42470 
.49135 
.56302 
.64088 
.72649 
.82214 
.93125 
.05932 
.21601 
.42049 
.71939 
.80006 
.89204 
.99889 
.12599 
28210 
.43272 
.75900 
.18758 
.05025 
.19820 
.36943 
.57151 
.81563 
,12031 
,51832 
07568 
95527 
75766 
06937 
43087 
85804 
37512 
02163 
86770 
05463 
93120 
78364 


DELC3 

.00000 
.01246 
.02286 
.03263 
.04205 
.05122 
.06022 
. 06909 
.07785 
.03653 
.09515 
.13771 
.13006 
.22281 
.26643 
.31133 
.35794 
.40675 
.45334 
.51347 
.57315 
.63381 
.71257 
.79782 
.90041 
.03189 
.22001 
.27006 
.32680 
.39229 
.46962 
.56384 
.68376 
.34700 
.09632 
.58629 
. 66395 
.76420 
.375 9 7 
.01021 
17659 
39211 
69076 
15535 
08588 
24428 
42714 
64212 10. 
90082 10. 
22206 11. 
63898 12. 
21780 13. 
12000 16. 
93054 20. 


D I FI 

.00000 
.00003 
.00014 
.00031 
.00055 
.00036 
.00125 
.00171 
.00225 
.00237 
.00357 
.00835 
.01544 
. 025 16 
.03737 
.05405 
.07425 
.09924 
.12996 
.16773 
.21433 
.27234 
.34566 
.44055 
.56798 
.74976 
.03899 
.12061 
21513 
.32653 
.46110 
62870 
.84706 
.15175 
.63035 
.60517 
.77323 
, 96800 
,19799 
,47619 
82380 
27841 
91581 
92313 
99100 
34902 
76436 
25530 
,84980 
59341 


OIF2 

DIF3 

PHI 

del: 5 

.OOCOO 

.OOOCG 

.00000 

.0G00Q 

-.00077 

. 0024o 

.01003 

.01CC0 

-.00157 

.00236 

.02014 

.02000 

-.00233 

.00263 

.03031 

.03000 

- .00306 

.30205 

.04055 

. ^4^Uu 

-.00374 

.00122 

.05086 

. u5 0 u d 

-.00436 

.00022 

.06125 

. 05000 

-.00493 

- .0009 4 

.07171 

. a f b 4 0 

- .Q05*o 

- .002*5 

.08225 

.02000 

-.00533 

-.00347 

.09237 

.09000 

- .00626 

- . 3C435 

.10357 

. : 0000 

- . 007'* 0 

-.01229 

.15835 

.15000 

-.00600 

-.01994 

.215^ 

.20000 

-.30265 

- . 0271 9 

.27516 

. 000 

.00329 

-.03357 

. J J F 8 / 

.30000 

.01224 

- . 03867 

. *0405 

.35000 

.02470 

-.04206 

.47425 

.^0000 

.04135 

-.04325 

.54924 

.*5000 

.06302 

-.0416o 

.62996 

.50000 

.09088 

-.03653 

.71773 

.55000 

.12649 

-.02685 

.31433 

.60000 

.17214 

-.01119 

.92234 

.65000 

.23125 

.01257 

1 . 34566 

.70000 

.30932 

.04782 

1.19055 

.75000 

.41601 

.10041 

1.36798 

.30000 

.57049 

.18189 

1.59976 

.35000 

.81939 

.32001 

1.93899 

.90000 

.89006 

.36006 

2.03061 

.91000 

.97204 

.40680 

2.13513 

.92000 

1 .06889 

.46229 

2.25658 

.93000 

1 . 18599 

.52962 

2.40110 

.94000 

1.33210 

.61384 

2.57870 

.95000 

1.52272 

.72376 

2.30706 

.96000 

1 . 78900 

.87700 

3.12175 

.97000 

1.20753 

1.11632 

3.61035 

.98000 

5.06025 

1 .59629 

4.59517 

.99000 

5.20720 

1 . 67795 

4.76423 

.99100 

5.37748 

1 . 77220 

4.96000 

.99200 

!. 57851 

1.88297 

5.19099 

.99300 

1.82163 

2.01621 

5.47019 

.99400 - 

.12531 

2.18159 

5.31880 

.99500 ■ 

■.52232 

2.39611 

6.27441 

.99600 ■ 

.07868 

2.69376 

6.91281 

.99700 ■ 

.95727 

3.15735 

7.92113 

.99800 - 

.75866 

4.08688 

9.99000 

.99900 - 

.07027 

4.24518 10.34812 

.99910 - 

.43167 

4.42794 10.76356 

.99920 - 

.85874 

4.64282 11.25460 

.99930 - 

.37572 

4.90142 11.84920 

.99940 - 

.02213 

5.22256 12.59291 

.99950 - 

.86310 

5.63938 13.56666 

.99960 - 

.05493 

6.21810 14.93353 

.99970 - 

.93140 

7.12020 17.09634 

.99980 - 

.78374 

8.93064 21.54219 

.99990 - 


d i f : : 

o.oe-oo 

:.3e-*o 

* . -c - '9 

? t:. '3 


DIFI2 

0. OE-OO 
6. ?E - "7 
5.6E- * 7 
-5 . 5 E- *6 
* 5 . 2E- ’5 
-3. 1c- 14 
•1-2E-13 
•4.3c- "3 
•1.3E- !2 


DIF:: 


0. OE-OO 


3.3E- 1 
-2.3E-1 
3.3E- 1 
a ’ 


0 . JC - U-* 

■“33 
03 


* i . . c - 
- 1 .6E- 
-2.2E- 
-3. 3E- 
-4.3E- 
•5 .3E- 
-6. 4 E- 


- 1 . 




-7.3c- 

-7.3E- 

-6.7E- 

-6.3E- 

-5 . 3E - 

-5 .2E - 

-4.5c- 

-3.5E- 

-2.5c- 

-1.25- 

-1. IE-03 

-9.3E-04 

-8.0E-04 

-6. 7E-G4 

-5.4E-04 

-4. IE-04 

■2. 95-04 

-1.75-04 

•6.9E-Q5 

■6.3E-05 

■5. IE-05 

4.3E-05 

■3.5E-05 

2.7E-05 

2.0E-05 

1.4E-05 

7.3E-06 

3.0E-06 


-2. IE-10 

-2. IE-09 

- 1 .2E-0S 

-5.3E-08 

- 1 . 3E - 07 

-5. IE-07 

-1.3E-06 

-2.9E-06 

-5.9E-06 

-1. IE-05 

-2.0E-05 

-3.2E-05 

-4.3E-05 

-6.6E-05 

-3.0E-05 

-7.7E-05 

-7.3E-05 

-6.7E-05 

-6.0E-05 

-5. IE-05 

-4. IE-05 

-3. IE-05 

-2.0E-05 

-9.5E-06 

-2.3E-06 

-1.3E-06 

-1.4E-06 

-1.0E-06 

-7.2E-07 

-4.7E-0 7 

-2.7E-07 

-1.3E-07 

-4.3E-08 

-7.9E-09 

-6.0E-09 

-4.4E-09 

-3.0E-09 

-2.0E-09 

-1.2E-09 

•6.7E-10 

■3.0E-10 

■ 9 . 9E - 1 1 

■1 .4E-11 


- 5 . 6E - " / 

-8.3E- -7 
-3.0E-15 
-3.5E- 4 4 
-2.9E- 13 
-1.9E-12 
-9.7E-12 
-4.2E- 11 
-1.5E-10 
-4.9E-10 
-1 .3E-09 
-3. IE-09 
-6. IE-09 
-9.2E-09 
-8.3E-09 
-8.0E-09 
-6.9E-09 
-5.5E-09 
-4. IE-09 
-2.7E-09 
-1.5E-09 
-6.2E-10 
-1.5E-10 
-8.7E-12 
-5.5E-12 
-3.2E-12 
-1.3E-12 
-8.7E-13 
-3.7E-13 
-1.2E-13 
-3.0E-14 
-3.9E-15 
-1. IE-16 
-1. IE-16 
0.QE+QQ 
0.0E+00 
0 . QE+00 
O.OE+OO 
0 . QE+00 
0.0E+00 
Q.OE+OO 
O.OE+OO 


0IFI4 

0 . OE-OO 
-+.2E- 17 
--+.5E- 17 
-3.3E- -3 
. 2E - ! 7 

- : . *E- 17 
-2.3E-17 
-2.3E - ■ 7 
-2.3E- 17 
-2.3E-17 
-5.6E- 17 

- 3 . oE - 1 7 
-5 .6E- 17 
-3.3E- 17 
-5 .oE- 17 
-1. IE-16 

0. OE-OO 

- 5 . 6E - 1 7 




17 

16 

16 

16 


-5.6E 
1 . IE 
-1. IE 
-1.1E 
-1. IE-16 
0. OE-OO 
0.0E-G0 
-1. IE-16 
-1. IE-16 
-1. IE-16 
-1 . IE-16 
0. OE-OO 
0. OE-OO 
0. OE-OO 
0. OE-OO 
0. OE-OO 
O.OE+OO 
O.OE+OO 
0.0E+00 
0. OE-OO 
O.OE+OO 
0.0E+00 
0. OE-OO 
0. OE-OO 
O.OE+OO 
0. OE-OO 
Q.OE+OO 
O.OE+OO 
0.0E+00 
O.OE+QO 
O.OE+OO 
O.OE+OO 
0. OE-OO 
0 . OE-OO 
0. OE-OO 
O.OE+OO 


0.2E-GO 
•5 . 6E- 1 7 
•5.0E-17 
0. OE-OO 
0.25-00 
0.25-00 
* 1 . ’ 5- *6 
0.25-00 
0 . „E-00 
0.25-00 
0.25-00 
0.25-00 
0.25-00 
0 .25-00 
0.05-00 
0.25-00 
0. 2E-00 
0. OE-OO 
0. OE-OO 
0.05-00 
0. OE-OO 
0.05-00 
0. OE-OO 
0. OE-OO 
0. OE-OO 
0. OE-OO 
0. OE-OO 
0. OE-OO 
0. 3E-00 
0. OE-OO 
0. OE-OO 
0. OE-OO 
0. OE-OO 
0. OE-OO 
0. OE-OO 
0. OE-OO 
G.OE+OO 
0. OE-OO 


EXACT CLELAND SOLUTIONS WITH DIFFERENCES BETWEEN CLELAND PHI 
SOLUTION SOLUTION APPLICABLE TO AND EXACT SOLUTIONS 

THIS PHI VALUE MARKED. 

NOTE THAT DEL VALUES OVER 
UNITY ARE UNPHYSICAL. 

PARAMETERS LEADING TO PHI VALUES 
OVER 1.6 CAUSE KORN2.EXE TO FAIL. 


5TH ERRORS FOR FIRST THRU 5TH ITERATIONS 
ITERATION NOTE THAT THE FIRST ITERATION IS VERY 

GOOD, AND THE THIRD IS EFFECTIVELY EXACT. 
THE CALCULATION USED DOUBLE PRECISION. 
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Page 


1C 


PROGRAM SOI VCU SC 
IMPLICIT REAL*8 (A-H ( 0-Z) 


FORMAT ( 1 COMPARISON OF CLSLAND "CURVE FITTING" SOLUTIONS 
‘FOR DELTA WITH EXACT 1 

2 \ x VALUE AND SIMPLE NEWTON ITERATION SOLUTION 1 // 

3 *14A9/) 


1 


WRITEvO, 10) 'DEL 1 , 'DEL Cl 1 , 1 DELC 
, 1 PH l 1 , 1 DEL 1 5 1 , 'DI F1 1 1 


2’ , ‘9ELC3* , 'D I FI 1 , ‘DIFZ 
, ' 0 1 F 1 2 • f • 0 1 F 1 3 1 , 'DIFI4 


1 C I F 3 1 
*0 1 FI5 1 


SElEC 


A RANGE OF EXACT DEL VALUES 3E7WEEN ZERO AND ONE 


DO IDEL=0,90 , 10 

DEL -I DEL /1 000. DO 
CALL SOLVE(DEL) 
ENDDO 


DO 1DEL=100,850,50 
DEL=IDEL/1 000 .00 
CALL SOLVE(DEL) 
ENDDO 


DO IDEL=900 , 990 , 10 
DEl=IDEL/1 000 . DO 
CALL SOLVE (DEL ) 
ENDOO 


DO IDEL=991 ,999, 1 
DEL=IDEL/100Q.D0 
CALL SOLVE(DEL) 
ENDDO 


DO IDEL=9991, 9999,1 
DEL=IDEL/ 10000 .DO 
CALL SOLVE(DEL) 
ENDDO 


END 


SUBROUTINE SOLVE(DEL) 

IMPLICIT REAL*8 <A-H,0-Z) 

DIMENSION D I FF(5 ) 

THIRD=1/3 .DO 

PHI=DEL/(1-DEL)**THIRD 

c THE PROBLEM IS TO SOLVE THIS EQUATION FOR DEL, WITH PHI KNOWN, 

C (WITHOUT USING THE KNOWN DEL VALUE, OF COURSE). 

C CLELAND "CURVE- FI TT I NG" SOLUTIONS APPLY TO DIFFERENT PHI REGIONS 
C AS INDICATED BY THE COMMENTS 

c if (PHI .le. .05) then 

DELC1 =PHI 

c else 

0ELC2= . 89086*PH I ** . 99301 
c if (DELC2.gt.0.6Q) then 

DELC3 = .6854*PH 1**0 . 8708 

c endi f 
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c end i f 

C NEKTON'S ITERATION SOLUTION 

C INITIAL APPROXIMATION IS EXACT SOLUTION Or THE EQUATION 
C WITH THE CUBES REPLACED 3Y SQUARES 

PHI 2=PH I *PH I 

DEL I =( SORT (PH 1 2*(4+PH 1 2) ) -PH 1 2 5/2 
C ITERATE 

C R IS THE FIRST RESIDUAL FORMULATION, AND RD IS ITS DERIVATIVE 
C USE S = 1-(1-R)**3, AND ITS DERIVATIVE SD, TO GIVE FASTER CONVERGENCE 
C AND AVOID PR08LEMS FOR LARGE PHI. 

DO 2 1=1,5 

RCOT=( 1 -DEL I )**TH IRD 

R=DEL I - PH I *ROQT 

RD=1+PH I*TH I RD*RCOT /< 1 -DEL I ) 

T= 1 - R 
T2=T*T 
S=1-T2*T 
SD=3*T2*RD 
DELI =0 EL I -S/SD 

C SAVE DIFFERENCE FROM EXACT SOLUTION, FOR DISPLAY 

2 DI FF( I )=DELI -DEL 

C URITE OUT THE SOLUTIONS AND ERRORS 

wn te(0 , 1 (9f9 . 5 , 1 P, 5E9 . 1 ) 1 )DEL , DELC1 , 0ELC2, DELC3 
1 , DELC1 -DEL, DELC2-DEL.DELC3-DEL, PHI , DELI ,DIFF 

end 
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APPENDIX C 

THE RCTREAC CCMPi -oR PRC GRAM 


The source programs RCTREAC . ?CR and CPICT.C 
following pages. The source programs, the executable 
and linked with Microsoft software) , and the commercial 



a single user license) have been supplied to . Rcoerc b .. --“1 

Kcmfeld of MSEC. The executable shculu net ce diS' — 

without a Graphic license. The executable RCTREACD.EXE - 

replacement dummy program CPLCTD.C, and dees not caou the -- --- — 
authorizes the general distribution of this cede. 


on 


(comp nee 
Graphic (with 
der and Dale 
to any user 
ode using the 
; package; RAI 


The program runs on IBM PC or AT compatibles, and can produce graphics 
on a very wide range of displays ard printers. Before runnrrg the program Lor 
the first time, it is necessary to run the program EQUIP. mo, *iuch comes with 
the Graphic software. This program allows the. user to select a cusp lay option 
and printer options from a menu of choices, and writes a permanent 
configuration file which is read each tire by the RCTFuAC program. More 
details are of course provided with the scrhwure. The environment . parameter 
GPC should be set to C:\GPC, or to seme ether DCS directory containing font 
files and the files written by EQUIP.EXE. 


The program solicits keyboard input for each . or eigr.t parameters 
defining the problem. For each parameter, the default is the value used the 
last time, which is saved in the file RCTREAC.DAT and used if the input data 
is null. Once each case is finished, the program begins again, soliciting new 
keyboard data. The sequence of cases is terminated using Cirl-Break (or any 
other interrupt such as invalid data, Ctri-C, power off, cr Aut— Ctrl-Cel ) at 
any stage, or by providing a negative time m response to the first prompt of 
a new case. 


For each case, the program writes its output to the screen. After 
producing its text output for each case, the program produces graphs. These 
are written both to the screen and to the graphics file RCTREAC. TKF. Printer 
output is produced by responding to the beep prompt with the keys L or 1 (for 
large full-page plots at high or low resolution), or M or m (for medium half- 
page plots at high or low resolution) . 

Once the run is over, the graphics file can be played back using the 

command 


PLAY RCTREAC. TKF 

(where PIAY.EXE is another Graphic program). This leads to . the same beep 
prompts as the graphs are displayed, with the same printer options. There is 
little point in using this option, however, since the RCTREAC program runs 
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quickly, and can be run again to obtain any desired graphical output. The 
full list of beep options is below. 


Sairtole screen output (cbrained using Ctrl-?) is also shewn below. 
File output is similar, except that the prompts and data entry are not shewn. 
Iters after the colons are entered from the keyboard. The scu^-s programs are 
listed on the following pages. 


The allowable choices are: 

c --> convert TKr file to anottier format 
l large, 1 ow “esolution plot 

L --> large, hign resolution plot 
m --> meai urn sized, low resolution plot 
M --> medium sized, high resolution plot 
p --> cnange printer selections 
q --> quit and close files 
v --> vary PostScript parameters 
w - - > shrink picture to fit window 
z --> zoom marked area to fill screen 
CR --> go on to next plot 
space bar --> to return to picture 


ENTER ANY CHANGES FROM THE BRACKETED DATA USED LAST. 

ALWAYS INCLUDE THE DECIMAL POINT WHEN YOU ENTER DATA. 

THEN PRESS THE ENTER KEY. 

TOTAL TIME IN HOURS (NEGATIVE TO QUIT) [ 10.0000] : 

PARTICLE DIAMETER (MICRONS) [200.0000] : 

ROTATION RATE (RPM) [ 2.0000] : 

REACTOR RADIUS (CM) [ 3.0000] : 

PARTICLE DENSITY [ 1.1000] : 

FLUID DENSITY [ 1.2000] : 

TEMPERATURE (DEG c)[ 20.0000]: 

ALTERNATE VISCOSITY (POISE) (NEGATIVE MEANS USE WATER § T) 


. 0100 ] : 


PROBLEM PARAMETERS: 
TIME PARTICLE 

(HR) DIAM (MICRON) 

10.0 200.00 


ROTATION RATE 
RPM RAD/SEC 
2.000 .2094 


REACTOR DENSITIES VISCOSITY 

RADIUS (CM) PARTICLE FLUID POISE 

3.0000 1.1000 1.2000 .01000 


SPIRAL ORBIT DESCRIPTION: 

FINAL AREA XOFFSET 

CONCENTRATION FRACTION (CM) 
2.01746 .21162 -1.0398 


BRADIUS S RADIUS 
(CM) (CM) 

1.9602 1.3801 


REYNOLDS TAYLOR 

number number 

2.6E-01 2.5E-03 


SPIRAL GROWTH 
PER RADIAN 
-4.7E-05 


ROTATION RATE CHOICE: 
OMAX OMIN PHI 

1/SEC 1/SEC 
.354 .072593 .43850 


DELTA 

.37493 


AREA 

FRACTION 

.21447 


ROTATION RATE 
1/SEC RPM 

.1936 1.8489 
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ROTREAC. 

FOR 


Friday, January 11, 1991 4:25 pm 

1 



INTERFACE TO SUBROUTINE 3GNPLOT [C , AL I AS : * _bgnD.ct ' ] 

Z 


X 

(MP , I SC , TKF I LE; 

3 



INTEGER *2 MP,ISC 

4 



CHARACTER *20 TKFILE [REFERENCE] 

5 



END 

6 

i 



INTERFACE TO SUBROUTINE CPLOT [C] < T I MH f AD I AM,Rr M , BRAD , DRHO, VI SCOS 

3 


1 

, ARE RAT , XSP I RAL , BRAD IUS , SRAD IUS 

9 


1 

f OMAX r CM I N , OMOPT , FCPT ) 

10 



REAL *4 TIMH,ADIAM,RPM, BRAD, DRHO, VISCCS 

11 


X 

f ARERAT , XSP I RAL, 3RAD IUS, SRAD IUS 

12 


X 

, CMAX , DM I N , OMOPT , FOPT 

13 



END 

14 




15 



INTERFACE TO SUBROUTINE STOPPLOT [C, ALIAS: ' jstccsiot '] () 

16 



END 

17 




18 



PROGRAM ROTREAC 

19 




20 

C 


ANALYZE ROTATING REACTOR PERFORMANCE 

21 

c 


ASSUME FLUID REMAINS IN SOLID-BODY ROTATION 

22 




23 



IMPLICIT REAL*4 (A-H f O-Z> 

24 




25 



LOGICAL LEXIST 

26 




27 



MP - 1 

23 



ISC = 1 g 1 

29 



CALL SGNPLOTCMP, ISC, 'ROTREAC. TKF'C) 

30 

31 

c 


GET INITIAL VALUES OF PARAMETERS FROM DATA FILE CR INITIAL DEFAULTS. 

32 

c 


OPEN DATA FILE. 

33 




34 



I NQUIREC FILE=' ROTREAC. DAT ‘ , EX IST=LEX 1ST) 

35 




36 



IF (.NOT. LEXIST) THEN 

37 



TIMH=20. 

38 



TEMP=20 . 

39 



RHOP=1 . 05 

40 



RHOF=1 . 00 

41 



VALT=- 1 . 

42 



AD I AM=20 

43 



BRAD=3 

44 



RPM=1 . 

45 



0PEN( 1 , FILE= 1 ROTREAC.DAT 1 , STATUS= 1 NEW 1 ) 

46 



ELSE 

47 



OPEN( 1 ,FILE=‘ ROTREAC.DAT 1 , STATUS= 'OLD 1 ) 

43 



READ(1 , 1 (/8F10.5 ) 1 ) 

49 


1 

T IMH, TEMP, RHOP,RHOF,VALT ( ADI AM, BRAD, RPM 

50 



ENDIF 

51 




52 

c 


READ IN CHANGES FROM DATA OF LAST RUN 

53 




54 

c 


FORMATS 

55 




56 

30 


FORMAT (1X,A,F8.4, ■] 

57 

31 


FORMATCF10.3) 

58 

32 


FORMAT ( / 1 ENTER ANY CHANGES FROM THE BRACKETED DATA USED LAST 

59 


1 

/' ALWAYS INCLUDE THE DECIMAL POINT WHEN YOU ENTER DATA. 1 

60 


1 

/• THEN PRESS THE ENTER KEY.') 

61 




62 

C 


CONTROL RETURNS HERE FROM THE END OF THE PROGRAM 

63 





Page 1 
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ROTREAC . FCR 


Friday, Januarv II, 1991 4:25 pm 


?3ge 2 


64 33 

65 

oo C 

67 

63 

69 

70 

71 

72 

73 

74 

75 

76 

77 C 

78 C 

79 

80 
81 
82 

83 

84 

85 C 

86 

87 

88 

89 

90 

91 

92 

93 

94 C 

95 

96 

97 

98 

99 

100 c 

101 c 

102 

103 

104 

105 

106 

107 

108 

109 

110 

111 c 

112 

113 

114 

115 

116 

117 C 

118 C 

119 

120 
121 
122 

123 

124 

125 

126 


WR I 75(0,32) 

READ TIME I N HOURS, STOP ON NEGATIVE. CONVERT "0 SECONDS 

WRI7E(0,30) ‘TOTAL TIME IN HOURS (NEGA* 7.E TO QUIT: [\TIMH 
READ (0,31) TIMHD 
IF (TIMHD .LT. 0.0) THEN 
CALL STOPPLOT 
STOP ' 1 

END IF 

IF (TIMHD .NE. 0.0) TIMH = 7 1 MuD 
TIM$cC=T;hh*3600. 

READ IN PARTICLE DIAMETER IN MICRONS 
COMPUTE THE PARTICLE RADIUS IN CM. 

UR I TE(0,30) ’PARTICLE DIAMETER (MICRONS) [' , AD I AM 
READ (0,31 ) AD I AMD 

IF (AD I AMD .GT. 0.0) AD I AM = AD I AMD 
ARAD=AD I AM*Q . 0001/2 . 

READ IN REACTOR ROTATION RATE IN RPM AND CONVEX ’0 RAD/SEC 

WR I TE(0 ,30) 'ROTATION RATE (RPM) [',RF* 

READ ( 0,31 ) RPMD 

IF (RPMD .NE. 0.0) RPM = RPMD 

ONE = 1 

PI=4*ATAN(0NE) 

OMEGA=RPM*2*PI /60 

READ IN REACTOR REACTOR RADIUS (CM). 

WRITE(0,30) 'REACTOR RADIUS (CM) [',BRA0 
READ(0 , 31 ) BRADD 

IF (BRADD .GT. 0.0) BRAD = BRADD 

READ IN SPECIFIC GRAVITY (DENSITY IN GM/CC) FOR ’HE FLUID 
ANO THE PARTICLES. 

WRITE(0,30) 'PARTICLE DENSITY [\RHOP 
READ (0,31 ) RHOPD 

IF (RHCPD .GT. 0.0) RHOP = RHCPD 

WR I TE(0,30) 'FLUID DENSITY [ ' , RHOF 
READ(0, 31 ) RHOFD 

IF (RHOFD .GT. 0.0) RHOF = RHOFD 

READ IN REACTOR WATER TEMPERATURE IN DEGREES CES’IGRADE (FOR VISCOSITY) 

WRITE( 0 ,30 ) 'TEMPERATURE (DEG c) [ 1 , TEMP 
READ(0 ,31 ) TEMPO 

IF (TEMPD .GT. 0.0) TEMP = TEMPD 

NATIONAL BUREAU OF STDS. FORMULAE FOR UATER VISCOSITY ( POISE, 

OR 1 GM/SEC CM); CRC HANDBOOK OF CHEM. AND PHYS., 70TH ED . , p. F-40. 

IF (TEMP. LE. 20.) THEN 

V=( 1301 . /(998. 333+8. 1855*(TEMP-Z3.)+0. 00585* 

((TEMP-20. )**2))> - 3.30233 
VISCOS=(10.0)**V 

ELSE 

V=( 1 .3272* (20 . - TEMP) - 0 . 001053*((’EMP-2Q. )**2) )/ 
(TEMP+105. ) 
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Pace 


127 VISCOS=. 01002 * (10.0)**V 

123 END I F 

129 

130 C PUT IN ALTERNATE VISCOSITY IN POISE IF DESIRED ,iN0 IGNORE TEMPERATURE; 

131 

132 WR I TE(Q ,30) 

133 1 1 ALTERNATE VISCOSITY (POISE) (NEGATIVE MEANS USE -~*ER 3 T) [' 

134 2 , VALT 

135 READ( 0,31 ) VALTD 

136 IF (VALTD .NE. 0.0) VALT = VALTD 

137 IF (VALT .GT. 0.0) 7ISC0S = VALT 

138 

139 C WRITE OUT DATA (DEFAULTS FOR NEXT RUN) 

140 

141 REWIND 1 

142 URITEO , 1 (8A10/8F10.5) 1 ) 

143 1 *TIMH 1 , ‘TEMP* # 'RHOP' , *RHOF 1 , ' VALT « , ' AO IAM • , 'BRAD 1 , l RPM‘ 

144 1 ,TIMH , TEMP # RHOP , RHOF , VALT , AD I AM , BRAD , RPM 

145 REWIND 1 

146 

147 C SPECIFY THE MEAN EARTH GRAVITATIONAL CONSTANT (9E0 3/CM**2). 

148 

149 C COMPUTE THE DRAG CONSTANT, CDRAG (1/SEC) WHICH IS THE RATE AT 

150 C WHICH MOTION DECAYS THROUGH DRAG FORCES ALONE. I' IS ONLY 

151 C USED MULTIPLIED BY RHOP; THE PRODUCT IS CALLED COR . 

152 

153 C COMPUTE THE OFFSET OF THE SPIRAL CENTER, XSPIRAL (CM) FROM 

154 C THE REACTOR CENTER, ALONG THE X AXIS. 

155 

156 C COMPUTE THE CENTRIFUGAL GROWTH RATE D (1/SEC) FCR THE SPIRAL RADIUS. 

157 

158 C THE FINAL RADIAL DISPLACEMENT FROM THE SPIRAL CENTER OF ANY 

159 C PARTICLE IS DELR TIMES THE INITIAL RADIAL D I SPLACEMENT . 

160 C THE PARTICLE CONCENTRATION STAYS UNIFORM, AND IS MULTIPLIED 

161 C BY CONC. 

162 

163 G=980. 

164 

165 CDR=9.*VI SCOS/ (2.*ARAD**2) 

166 CDRAG=CDR/RHOP 

167 

168 DRHO= (RHOP- RHOF ) 

169 ADELRH=ABS(DRHO) 

170 

171 XSPIRAL=(G*DRHO)/(OMEGA*CDR) 

172 AXSP=ABS(XSPI RAL ) 

173 

174 D=(DRHO*QMEGA**2)/CDR 

175 

176 ETA = DMIMSEC 

177 AETA = ABS(ETA) 

178 DELR=EXP(ETA) 

179 CONC=1 / (DELR*DELR) 

180 

181 C COMPUTE THE REYNOLDS AND TAYLOR DIMENSIONLESS PARAMETERS. 

182 C BOTH SHOULD BE MUCH LESS THAN UNITY FOR A VALID SOLUTION. 

183 C THE RATIO D/OMEGA (SPIRAL GROWTH PER RADIAN) MUST ALSO BE 

184 C SMALL FOR A VALID SOLUTION. 

185 

186 REY = RHQF*ARAD*AXSP*QMEGA/V1SCQS 

187 TAYLOR = RHOF*ARAD**2*OMEGA/VI SCOS 

188 SPGRPR = D/OMEGA 

189 
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Page 


190 

C 

191 


192 

c 

193 

c 

194 


195 

c 

196 

c 

197 


198 


199 


200 


201 


202 


203 

c 

204 


205 

c 

206 


207 

c 

208 

c 

209 

c 

210 


211 


212 


213 


214 

c 

215 

c 

216 

c 

217 


218 


219 


220 


221 


222 


223 

c 

224 

c 

225 


226 

c 

227 


228 

c 

229 

c 

230 

c 

231 

c 

232 


233 


234 


235 


236 

c 

237 

c 

238 

c 

239 

c 

240 


241 


242 


243 


244 


245 


246 


247 


248 


249 

2 

250 


251 



252 


COMPUTE THE RAOIUS RSP'.RAL (cm) OF THE SMALL CIRCLE. 

HEAVY PARTICLES ORIGINALLY IN THE CIRCLE DO HOT SPIRAL OUT TO 
HIT THE WALL IN TIME TIMH. 

BUOYANT PARTICLES WHICH 00 NOT HIT THE WALL DURING^HE FIRST 
SPIRAL SPIRAL INUARDS TO INSIDE THIS CIRCLE IN ".ME TIMH. 


DELTA = AXSP/BRAO 

BRAD I US = 3RAD-AXSP 

SRADIUS = BRAD IUS’EXPL -AETA) 

ARERAT = ( SR AD I US/BRAD )**2 

THE REMAINING CALCULATIONS IGNORE THE GIVEN ROTATION RATE 

COMPUTE THE MAXIMUM AND MINIMUM VALUES OF THE ROTATION RATE (RAO/SEO. 

THE MINIMUM IS DEFINED SO THAT THE SPIRAL CENTER IS AT THE REACiOR RADIUS 
THE MAXIMUM IS DEFINED ARBITRARILY BY ABS(ETA)=1, SO THAT PARTiCLES 
MIGRATE INUARDS OR OUTUARDS CENTR I FUGALLY 3Y A "ACTOR e IN TIME TIMH. 


OMAX = (CDR/(ADELRH»TIMSEC))**0.5 
OMIN = G’(ADELRH)/(BRAD*CDR) 


AN OPTIMIZED ROTATION RATE IS COMPUTED 3Y MAXIMIZING THE RATIO OF THE 
AREA INSIDE THE SPIRAL RADIUS TO THE REACTOR AREA. in>nllulln „. 

FIRST DEFINE PHI. NOTE THAT PHI**3/2 = (OMM I N/0MMAX)**2 = (ARAD/AMAX) 6 


THIRD=1/3.D0 

PHI = (2.*ARAD**6*G**2*TIMSEC/BRAD**2)**THIRD 
*(2.*A0ELRH/(9.*VISC0S) ) 

SOLVE DEL = PHI*(1-0EL)**(1/3) 

DEL IS OMOPT/OMMIN = XSP I RAL/BRAD 


NEWTON'S ITERATION SOLUTION 

INITIAL APPROXIMATION IS EXACT SOLUTION OF THE ECUATION 
UITH THE CUBE ROOT REPLACED BY A SQUARE ROOT. THE POSITIVE 
ROOT IS CHOSEN, AND FORMULATED TO AVOID ROUNDING ERROR FOR 
VERY LARGE PHI 

PHI2=PHI*PHI 

DELI=PHI2*2/(SQRT(PHI2*(4+PHI2))+PHI2) 


R IS THE FIRST RESIDUAL FORMULATION, AND RD IS ITS DERIVATIVE 

USE S = 1-(1-R)**3, AND ITS DERIVATIVE SD, TO GIVE FASTER CONVERGENCE 

AND AVOID PROBLEMS FOR LARGE PHI. 

DO 2 1=1,5 

ROOT =(1 -DELI )**TH IRD 
R=OELI -PHI*ROOT 
RD=1+PHt*R00T*THIRD/( 1 -DELI ) 

Q=1 -R 
Q2=Q*Q 
S=1 -Q2*Q 
SD=3*Q2*RD 
DELI=DELI -S/SD 

DELOPT=OELI 

FOPT = ( 1 -DELOPT)**2/EXP(PHI**3/DELOPT**2) 
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Page 5 


253 

254 

255 

256 C 

257 C 

258 

259 

260 
261 
262 

263 

264 

265 

266 
267 
263 

269 

270 

271 

272 

273 80 

274 

275 

276 

277 

278 

279 

280 81 
281 
282 

283 

284 

285 

286 

287 83 

288 

289 

290 

291 

292 

293 

294 82 

295 

296 

297 

298 

299 

300 

301 

302 

303 

304 

305 

306 

307 

308 

309 

310 

311 C 

312 

313 


CMOPT = (G/(2.*SRAD’TIMScC*( 1 -DELCPT)))**"HIRD 

AN APPROXIMATE UPPER LIMIT ON PARTICLE RADIUS, AMAX (CM), IS COMPUTED 
3Y EQUATING OMAX AND OMIN. ALL OTHER CALCULATIONS JSE THE GiVEN SiZE. 

AMAX = (9.*VISCOS/(2.*ADELRH))”0.5* 

1 (BRAD”2/(G*»2*T;MScC))**C1./6.) 

OMAX = AMAX*20000. 

UR I TE(0, SO )T I MH, AD I AM, RPM, CMECA, 3RAD, RrC=, RHCF, VI SCCS 

IF (BRAD.GT.ABS(XSPIRAL)) THEN _ 

UR I TEC 0,31 )CONC , ARERAT ,XSP I RAL , 3R AO I US, SR AD IUS, REY , TAYLOR , i.-jKrR 

else 

WRITE(0,83)C0NC, ARERAT, XS?IRAL,3RADIUS,SRADIUS, REY, TAYLCR,S?GRPR 

END IF 

UR I TE (0,32) OMAX , OMIN , PH I ,DELCPT, FOPT,CMC?',OMOPT*30/PI 

FORMATC// 1 PROBLEM PARAMETERS: 1 ,/, 

1 1 TIME PARTICLE ROTATION RATE 

1 1 REACTOR DENSITIES VISCOSITY 1 ,/, 

2 1 (HR) DIAM(MICRON) RPM RAD/SEC ', 

2 1 RAD I US ( CM) PARTICLE FLUID POISE 1 ,/, 

3 F6.1,F12.2,F10.3,F3.4,F12.4,F9.4,F8.4,F9.5,//) 

FORMATC SPIRAL ORBIT DESCRIPTION: 1 ,/, 

1 > FINAL AREA XOFFSET BRAD IUS SRADIUS ', 

1 'REYNOLDS TAYLOR SPIRAL GROUTH 1 ,/, 

2 1 CONCENTRATION FRACTION (CM) (CM) (CM) 1 , 

3 1 NUMBER NUMBER PER RADIAN 1 ,/, 

4 F10.5,F12.5,F8.4,F9.4,F8.4,1P,2E9.1,E11.1//) 

FORMATC SPIRAL ORBIT DESCRIPTION IMPOSSIBLE: 1 ,/, 

1 ' FINAL AREA XOFFSET BRAD IUS SRADIUS 

1 'REYNOLDS TAYLOR SPIRAL GROUTH 1 ,/, 

2 1 CONCENTRATION FRACTION (CM) (CM) (CM) ', 

3 1 NUMBER NUMBER PER RADIAN 1 ,/, 

4 F10.5,F12.5,F8.4,F9.4,F8.4, IP, 2E9. 1 ,E1 1 . 1//) 

FORMATC ROTATION RATE CHOICE: 1 ,/, 

1 1 OMAX OMIN PHI DELTA AREA ', 

1 'ROTATION RATE ',/, 

3 1 1/SEC 1/SEC FRACTION 1 , 

1 1 1/SEC RPM 1 ,/ , 

5 F7.3,F9.6,3F9.5,2F9.4,//) 

URITE (0,*) 'ENTER N IF YOU DO NOT UANT A PLOT 1 
READ (0, '(AD 1 ) CNO 
IF (CNO.EQ. 1 N 1 ) GO TO 33 

CALL CPLOT (TIMH,ADIAM, RPM, BRAD, DRHO.VISCOS 
1 , ARERAT , XSP I RAL , BRAD I US , SR AD I US 
1 ,OMAX,OMIN,OMOPT , FOPT) 

GO TO 33 

IT KEEPS SOLICITING A NEU CASE, AND STOPS ON A NEGATIVE TIMH VALUE 
ENO 


!§ 


OF PQO% QUALITY 
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CPLOT. C Friday, January 11, 1991 4:24 pm Page 1 

1 ^include "graph i c. h" /* Will include all necessary include files */ 

2 

3 /* 

4 CALL CPLOT (TIMH , AD I AM , RPM , 3RAD ,DRH0, VI SCCS , ARERAT 

5 1 , XSP I RAL , BRAD I US , SRAD IUS 

6 1 , OMAX , 3M I N , OMOPT , FOPT ) 

7 */ 

3 

9 void cplot (float timh, float adiam, float rpm, float braa, float arho, float viscos. 


10 


float arerat, float xspi 

ral, float bradius, float sradius, 

11 


float omax, float omin, 

float omopt, float fop: ) 

12 

13 

< 



14 

/* 



15 
1 >4. 


declarations 


j D 

17 

* / 

float brads, xs, bra, sra,rgrowth; 


18 


int i ,nxdiv,nydiv,npts; 


19 


char buffer [993; 


20 


float omC301], y[301], z [301 j # w C301] ; 


21 


float ymin, ymax, ystep, omaxd; 


22 




23 

r 



24 


first plot 


25 

V 



26 


if (abs(xspi ral )<brad) 


27 


< 


28 


/* CONVERT TO RADIUS 2.3 INCHES */ 


29 


xs=(2.3/brad)*xspi ral; 


30 


bra=(2.3/brad)*bradius; 


31 


sra=(2.3/brad)*sradius; 


32 


brads=2.3; 


33 


startplot(7); /* initializes each plot and sets the background color */ 

34 


/* 0-7 are black, blue, green, cyan; 

red, magenta, brown, white */ 

35 


font (3, "simp lex. fnt", 1 \3 1 0 1 f "triplex. fnt 

", 1 \3 1 1 ’ , "comp l ex. fnt", 1 \3 1 2 ' 

36 


ii M i i\. 
t i 


37 



/* loads your chosen fonts */ 

38 


rotate<0); 


39 


page(9. ,6.855); 

/* sets the output page size */ 

40 


setscale(O); 


41 


/* this is the same aspect ratio as 4095 by 3119 */ 

42 


/* cross(l); 

crossed axes if either axis goes through zero */ 

43 


/* area2d(8.0,5 .5); sets the 

area of the plot (drawing page) */ 

44 


color(2); 

/* Green */ 

45 


box ( ) ; 


46 


col or (5 ); 

/* sets color to magenta */ 

47 


tmargin( .15); 


48 


ct l i ne("\31 1 ROTARY REACTOR CROSS SECTION 

", .31 ); /* titles plot */ 

49 


color(4); 

/* sets color to red */ 

50 


dashf ( 1 ); 


51 


c i r c l e( 2 . 7+0 .,2.7+0., brads , 0 . , 6 . 3 ) ; 


52 


dashf (9); 


53 


c i re l e( 2 . 7+xs , 2 . 7+0 . , bra , 0 . , 6 . 3 ) ; 


54 


ci re le( 2. 7+xs, 2. 7+0. ,sra,0.,6.3); 


55 


color(2); 

/* Green */ 

56 


datei t( * \3 1 0 1 ); 


57 




58 


page(3. 2,4.6); 


59 


pgshi ft(5.4, .4); 


60 


l inesp( 1.8); 


61 


setscale(l); 


62 


box ( ) ; 


63 
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Lmargin( .2); 

tmargin( .0) ; 

ctl ine("\31 1DATA", .7); 

tmarginC 1.5); 

sprintf (buffer, "%S. If", timh ); 

r t L ine(buffer, .45); 

s pr i n t f ( bu f f e r , "%3 . 1 f M , aa i am ) ; 

rt L ine(buffer, .45); 

spr i ntf (buff er , "7=8. 1 f ", rpm) ; 

rt l i neCbuf fer, .45); 

spri ntf (buf fer, "%8. 1 f " , brad) ; 

rtl ine( buffer, .45 ); 

spri ntf ( buf f er, "%8.3f " , drho) ; 

rtl ine< buffer , .45); 

spri ntf (buf fer, "%8. 4f ", vi scos) ; 

rtl ine(buf fer, .45); 

tmargi n(1 .5); 

It l ine( "Hours", .45); 

It l ine( "Mi crons", .45); 

Itl ine("RPM", .45); 

It l ine("Reactor Radius" , .45 ) ; 

It l ine( "Densi ty Excess", .45); 

Lt l ine< "Viscosity", .45); 

tmargi n(7.0) ; 

ct L i ne("\31 1 RESULTS" , . 7) ; 

tmargi n(8.5); 

spri ntf (buf fer , H %8.4f",arerat ) ; 
rt l i ne(buf fer, .45); 
spr i ntf (buf fer , ’'%8 ,4f " , xspi ral ); 
rt l ine(buffer, .45); 

if (drho>0.) rgrowth=bradius/sradius; 
else rgrowth=sradius/bradius; 

sprintf (buf fer, "%a.4f",rgrowth) ; 
rt l ine(buf fer, .45) ; 

tmargin(8.5); 

Itl ine("Area Fraction", .45); 

Itl ineO'Spi ral Offset", .45); 

Itl ine("Radius Growth", .45) ; 

endplotO; 

) 

/* END OF IF BLOCK ON XSPIRAL TEST. START OPTIMIZATION PLOT */ 

startplot(7); /* initializes each plot and sets the background color V 
/* 0-7 are black, blue, green, cyan; red, magenta, brown, white */ 
f on t ( 3, "simplex, fnt", - \31 0 • , "triplex, fnt", '\311 "complex. fnt", *\312 ' 

llll I I J . 

/* loacs your chosen fonts */ 

rotate(90); 

page(6.S55,9. ); /* sets the output page size / 

setscale(O); 

/* this is the same aspect ratio as 4095 by 3119 */ 
color(2); /* Green V 

box(); 

color(5 ); /* sets color to magenta V 

tmargi n( .12); 

ct l i ne("\31 1 ROTATION RATE OPTIMIZATION", .25); /* titles plot / 
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color(4); .* sets color to red V 

datei t ( 1 \310 1 ) ; 


ymi n=0. ; 
ymax=1 ; 


ystep= .2; 
page(6.355 ,5.5); 
pgshift(0. ,2.7); 
area2d(5. 355,4.0); 
coior( 10) ; 
grid( 1 ); 
fgrid(2, 2, 2) ; 
framed , 1 ); 
cross(0) ; 
fntchgC 1 \31 0 * ) ; 
omaxd=max(2*omax , 5*omi n) ; 
for (i=0; i <=300; i++H 


/* sets the outcut page size */ 

/* Plot diminisions */ 
/* 3lack grid */ 
/* Oottea l ine grid •/ 

I* SuDdivide the c.ct with tick marks only *7 
/* Draw a frame around the plot */ 
/* Do not cross axes */ 
/* Charge fonts for axes labels V 

/* Make vectors to plot */ 


om[i ] =oniin*sxp( log(omaxd/omin)*i/300. ); 
y [i] =1 -omin/om(i] ; 

z [i]=exp( * (om[i)/omax)*(om[i] /omax) ); 

wCi]=yCi]*zCil; 


> 

npts=301 ; 

xlog(omin,omaxd, 1 .1 f " , ymi n,ystep, ymax) ; 

colord 1 ); 

xname("\310RPM" ) ; 

yname("\310Radius Ratios"); 

color(IO); 

curve(om, y, npts, 0); 
color<14) ; 

curve(om, z, npts, 0); 
color (1 2) ; 

curve(om, w, npts, 0); 


/* Size of the x i y vectors */ 

/* Manual scaLed om-axis */ 

/* om-axis name */ 
/* Y-axis name */ 

/* Draw curve with no symbols */ 

/* Draw second curve */ 

/+ Draw third curve *f 


page(3. 2,2.2); 
pgshi ft( . 15, . 15); 
l inespC 1 .8); 
box(); 


ImarginC .2); 
tmargin( .0); 
ct l ine<"\31 1DATA" , .20); 


tmargin( .38); 

spr intf (buffer, "%8. If", timh); 

rt l ine(buf fer, .13); 

spr intf (buffer , "%8.1 f ", adiam) ; 

rt l ine(buffer , .13); 

spr intf (buffer, "%8. If", brad); 

rt l ineCbuf fer, .13); 

spr intf (buf fer, "%8.3f ",drho) ; 

rt l ine( buffer, .13); 

spr intf (buf fer, "%8.4f ", viscos); 

rt l ine(buf fer, .13); 


tmargin( .38); 

It l ineC’Hours", . 13); 

Itl ine< "Mi crons", .13); 

It l ineC'Reactor Radius", .13); 

It l ine< "Densi ty Excess", . 13) ; 

Itl ine<"Viscosi ty", . 13); 

page<3. 2,2.2); 
pgshift(3.5, .15); 

Ot 


w fS 

QUALITY 


- 30 - 


Rotating Reactor Studies 


Roberts Associates 


CPIOT.C 


Friday, January 11, 1991 4:24 pm 


Page 


190 

191 

192 

193 

194 

195 

196 

197 

198 

199 

200 
201 
202 

203 

204 

205 

206 

207 

208 

209 

210 
211 
212 

213 

214 > 


box ( ) ; 


lmargin< .2); 
tmargin( .0) ; 

ctline("\31 1 RESULTS”, .20); 


tmargin( .38); 

spr i nt f ( buffer , "%8 .3f omi n) ; 
rt l ine(buf fer, . 13) 
sprintf (buffer, "%8 
rtl ine(buffer, .13) 
sprintf (buffer, '*%8 
rt l ine(buffer, .13) 
sprintf (buffer, M %3 
rt L ine(buf fer, .13) 


3f " , omax ) ; 

3f " , omopt ) ; 

3f ", sqrt( f opt ) ) ; 


tmargin( .38); 

l tl ine ("Mini mum RPM",.13); 

Itl ine( "Maxi mum RPM»',.13); 

It L ineC'Optimum RPM M ,.13); 

1 1 L ine( "Optimum Fract ion'*, .13); 

endpLot( ) ; 

scrt(); /* text mode for interactive stuff */ 


CPLOTD.C 


Friday, January 11, 199 1 ^ ; 25 pm 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 


„ id mow 

float omax. float omin, float omopt, float fopt ) 


void bgnplot (int mp, int isc, char * tkfile) 

C 

> 

void stopplot () 

C 

> 


Page 1 
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PARTICLE orrtts IN A ROTATING LIQUID 


An earlier version of the appended report, with William Fcwlis as the 
first author, was submitted in early 1988 to Prof. Owen Phillips at Johns 
Hopkins University, for publication in the Journal of Fluid Mechanics.. A 
favorable report was received from the two referees, with only minor 
corrections required. Submission of a corrected version was decayed by the 
increasing incapacity of Dr. Fcwlis, who went on extended sick ^ , e 

summer of 1988, later took medical disability retirement, and died at the 

beginning of 1989. 


With the concurrence of Dale Kcmfeld and of Dr. Robert S. Snyder 
(supervisor of Dr. Fowlis at MSEC until his retirement) Dr. Roberts of RAI 
agreed to become first author and to resubmit the paper. The new version 
follows. 
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Roberts Associates, Incorporated 
380 West; Maple Avenue, Suita L-1A 
Vienna, VA 22180-5616 

(703)242-2115 cr (703)356-6520 
(Washington Metropolitan Area) 


Prof. O.M. Fhillips 

Journal of Fluid Mechanics 

Dept. of Earth and Planetary Sciences 

The Johns Hopkins University 

Baltimore, MD 21213 

Dear Professor Phillips: 

I enclose a revised manuscript; ror cur paper "Parcac-e Orbits m a 
Rotating Liquid". I believe this revision satisfies all the comments both of 
referees A and B. 

Bill Fcwlis died in December, 19S3 . He had been increasingly sick ^cr 
three years, with what was ultimacely diagnosed as a form of Parkinson's 
disease. He stepped wc riling in the summer of 1983. 

Dale Komfeld and Dr. Robert S. Snyder (of Code ES76, MSFC, ^Teh 
205/544-7813 and 7805) have agreed that I should be identified as first 
author, as reflected on the amended title page. Dr. Snyder is the contracoing 
officer's technic al representative on my company's contracts, and is 
Komfeld 's supervisor. He was also Bill Fcwlis 's supervisor urxil his death. 

The original manuscript was lost from the word-processor. You will 
note that this manuscript is a photocopy, with the changes cut-and-pasted in. 
I hope this is acceptable; it was easier than retyping the whole thing would 
have been. I do not have originals for the figures 1 and 2; so if you do mot 
have them, please let me knew and we will get replacements to you. 

I believe t his revision satisfies all the comments and suggestions of 
referees A and B. 

We have incl uded the references mentioned by referee A, and one ether 
new reference. Also, we have corrected the Ekman number definition; we thank 
him for pointing out this error. 


With regard to the first paragraph of referee B, our process has been 
used to make large spherical latex microspheres with quality and uniformity 
approaching that of microgravity processing, as described in Section 1, at the 
bottom of page 3. And our analysis has been used in setting the. rotation rate 
range for new rotating apparatus designs; we have expanded Section 5 sligntly 
to clarify this. 


We have modified the paper in response to suggestions 2 and 3 of 

referee B. With regard to suggestion 1, equations (27) through (29) are 
derived from equation (16) , using the approximation (24) , as sooted at the 
bottom of page 11. 

SuGoesoicn 4 is answered by the wr.cie of Section 4, where a spe^^ic 
optimization problem is solved. We have mcditied the text on page^io to 
emchasizs than this is dene without using (33). A siigbcr/ ^ di_*t=rent 
optimization problem michc have a slightly different answer for the optimum 
rotation rata. There is no strict upper limit on the rotation rata, bur the 
fas tar the chamber spins, the mere rapidly trie particles cantriruge outward -~o 
the walls or inward towards the axis. 

We have added Ficure 3, to clarity the optimization prcclem solved in 
Section 4 . And we have added brief ramarxs on page 25 concerning secondary 
flews driven by the excess density distribution due to variations u. the 
particle concentration. 

We thank the referees for their work. 

I apologize for the delay in retiming this paper. This should have 
been done much more quickly after Bill's death. 

Yours sincerely, 


Dr. Glyn 0. Roberts 
President 
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2. FORMULATION AND SOLUTION OF THE PROBLEM FOR PARTICLE ORBITS 


Basic mechanisms for the suspension of hea 
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.icles in a 
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rotating cylinder of fluid were discussed by Otco anc Loren 
but a formal solution was not obtained. A more complete study of 
the oarticle orbit problem was presented by Schacc, / but tne 
centrifugal buoyancy was not handled correctly. Related s^udxe- bv 
Dill and Brenner 8 , by Nadim, Cox and Brenner 9 , by Aoki, Shirane, 
Tokimoto and Nakagawi 10 , and by Annamalai and Cole , addressed pa- “ s 
of the problem. In their later work, Annamalai and Cole obtained 
an orbit solution essentially equivalent to ours, in the context 
bubbles. None of these references analyzed the optimization problem 

studied here. 

In this section the 

particle orbit problem is correctly solved. The ambient fluid is 
taken to be in solid-body rotation about the horizontal z axis 
through the origin (see Figure 2) with rotation rate n. Then the 
fluid flow in the (x,y) plane is given by 


u = n x x = (-Qy, ax) . 


( 1 ) 


The corresponding pressure distribution is 

d = d + p f (i a 2 r 2 - gy) , (2) 

where p Q is a constant, p - is tne fluid density and r is 
x 2 + y 2 , so that r is the cylindrical radius coordinate. The 
gradient of this pressure balances the gravitational .and 

centrifugal forces on the fluid. 

The vector equation of motion for a spherical particle of 
radius a, density p f volume V, and mass M at position _>L is 
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The g term is tne familia- ...wfl-— - . 

tea is the corresponding inward press-are gradient which would 

cancel the centrifugal acceleration cf the liquid which the 

. . . . was aooarentlv omitted by 

particle has displaced. -his t_m was a.. 

7 

Schatz . 

a the Stokes slow— flow 
The drag force is written us.ng 

approximation. In this viscous limit 


D = - 6-ta(x - 8 x x) 


( 


-r^p d^aa force opposes the 
where ti is the fluid visccsiuy. • - 

„ . , . < ^ fluid Th^ aDDroxirciat ion 

motion of the particle relative to tne flu . --- -- 

, , =ri( j *-he Taylor number be much 

requires that the Reynolds number and the / 

smaller than unity, where 


7 


VO 



Re = p-aV/n 


( 7 ) 


(3) 


Ta = P - a ^ ' 


. t ■ _ _ - ^ ~ ^ ' -O (~\J — | X i» X X| ) . 

and V is the relative motion uj. tn- -> -- -- I— — 

• ‘■■jp -.on-t'-iviai cornooner.ts 

Combining these expressions, — 

of the equation of motion can be written as 


" . Pf 2 


x + cx + — a x + cftv - 0 , 

°o 


n i P f 2 

y+cy+— ay- cnx = 


A Q 


(9) 


( 10 ) 


where the drag constant. 


2 o a 
P 


2 r 


(ID 


is the rate at which motion decays through drag forces alone and 
ip = P - p f . The term on the right hand side of equation (10) 
arises P from the particle weight less the Archimedes buoyancy 

force. The S! 2 terms on the left hand sides of the equations (9) 

_ • y- nr*o c c’ i in equation 

and (10) arise from tne centri-ugd or- 

(5). 

The coupled x and y equations can be solved conveniently by 
introducing the complex variable 


w = x + iy. 


( 12 ) 


8 



Multiplying equation (10) by i and adding equat: 
come lex, second-order, linear equ = -~on 


w + cv 


p f 9 

(— 




A o 


•Th6 Gsners- sc 

;icular intecra. 


:n cf equation (13) is 
a ccrolerentary funct: 


w = w 


+ Ai exp^t) + - 1 - 2 ex?U 2 t). 


For this problem the p< 


: ^ t ^ color integral is tne co: 


w o = x o 


A OQ 

W Q ” p cH 

D 


n __ _ 

. H f a.-i 

' 


The real and imaginary pa: 


•ts reoresent the equili 


the x-y plane. The complex con 
rr 0 are the roots of the quadratic auxiliary e«j. 


of the particle in 
‘‘2 


m 


^ + cm + ( — — ft -left) 0 , 
P. 


The complex constants and A 2 are determined 
values w ( 0 ) and w(0) at time zero, using the equa: 


A^ + A 2 — w( 0 ) — Wq , 

i 

m-j_A]_ + m 2 A 2 = w(0). 


(9) gives the 

( 12 ) 

rim. of a 

(14) 

- s t a n t 

( 15 ) 

rium position 
stants m 2 and 
on 

(16) 

the initial 
.ions 

(17) 

(13) 
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These simultaneous ecu a t 


?he general solution 


ions can be easily solved for Aj and Aj 

n (14) represents the superposition of 


;wo SO 


iral motions about t 


the equilibrium dcs-l-i 


'n (x Q , y Q ). 


Consiae: 


the one solution 


w = x +iy = 


An exo (m- t ) , 

1 “ X 


= aexp(i9) , 


m-]_ = m r + imj_ 


( 20 ) 


(21) 


Here a 


is the positive amplitude of A 1 and 


9 is its phase r while 


m r and rr^ are the real and imag 
imaginary parts of w are 


m t 

= a e r cos(m.t + 0) 


inarv part of m^ 


. The real and 


( 22 ) 


y = a e 


m t 

: r sin(rmt + e) 


(23) 


m t 

The instantaneous spiral radius is a • - and grows or decays 

exponentially depending on the sign of The phase angle ^t 

+ , determines the direction of the vector displacement (*,Y> » 

the x-y plane . 

We now use the inequality 


10 



( 24 ) 

n < < c . 



w. 


= x. 


1 V = 

- o 


C A 3 
ftC 3 


2ca 2 A: 
9 nft 


(25) 


so that the spirals are cer. 
disolaced horizontally z r on 
the drag force balances tne 
must be inside the rotating 
outside or on the rotating 
will strike the wall. This 


tered about 

a point on 

the 

x-axi 

3 , 

the axis or 

rotation . 

At 

this 

me i .o t 

* 


n. 0 1 - weicht . Natura.iv , tnis point 
reactor, since if tne point is 
reactor wall all the particle orbits 
outs a lower limit on C ; 


a 


> n 


min 




(26) 


where b is the radius of the rotating reactor. For 
| A p|/p ~ 0.1, g ~ 10 3 cm/sec 2 , b - 3.0 cm and the above value 

for c, a . is 7 . 4x10” 4 sec -1 , cf . Table 2. 
mm 

Using the assumption (24), the two roots of the quadratic 
equation (16) can be approximated as 


m-j_ = ift + d , 


(27) 


re 


Or KrfW 0UAU7Y 
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( 2 ~), 


= -C “ iJi r 


he growth or ceoey r-. 


, . £4. — *. <--n -••>-= i , in equation 


for tne nr 


Ad n~ 

0 = h “ ' 


which is real and muc: 


than n. Net 


that for the s; 


miral, the real part is c, so 


that the radius decreases 


exoo 


nentially on the very s 


hort drag tine 


scale, 1/c . Thus the 


, . 1 , . • u 0 ~~~-os ne'-licibly small in a few drag 

second soirsl solu< ^ 


ime scales, for any realistic mit-al con. 


itior.s. The ful. 


-eneral soluti ion tnars . o - a 


:an be 


■oxidated as 


x Q 4- ae 


cos ( 3t 


e) , 


(30) 


= ae dt sin ( at 4 e) , 


where a and a are determined by the initial displacement from the 
center of the spiral - 

Note that the radius of the displaement from the spiral 
center is «e dt , where d is much less than a and has the sign 
of ip. The rotation rate of the spiral is a, so tne particle 
rotates about the center of the spiral at the same rate as the 
fluid rotates about the axis. The exponent for S r°“th or decay 
of the radius during the experimental time T is 
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Whether Ap is p 


■'ns exccnen: 


:st of the pa: 


or they spira_ 


we recuire 


* <-* -h 

A 2 I .1 i 


; ve or n 


cr:er 


"’73 r W — 


\S sriral cur 


» f- ~ 1 lO ’ 7 C — ..wC!I< - — - — — v 


zf the reactor 


i 33 i ve . - nus 


2 2 

n < a . 


2a T : A - 


For T = 10 3 sec, or slightly larger tn 


’ j the u - a 6 r 


Physically, the spiralling in or out is caused by the 


:en tr i fucai 


' ^ 3 s s u ra 


, . r , This force is oalancec c y tee 

That is why it has the sign - 3 * 


inwards cr 


outwards 


By equating the lower and upper limits, ( -o ) 
for a, we can obtain an approximate upper limit or tne r— 
the particles for the rotating reactor to be userul, 


6 _ r 9n ^3 Jl 

lax ^ 2 I Ao t 2 . 


( 34 ) 


Using n ~ 10 -2 gm/cm sec, |ad| =0.1 gm/cc, b - 3 cm, g 
cm/sec 2 and T = 10 5 sec, as before, we obtain a_ ax as 1.— * 10 
cm, corresponding to a diameter of 284 pm. ihis upper 1-u.t 
further discussed in the following section. 

• - .-,r 
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4. OPTIMIZATION OF THE ROTATION 


From the previous section, tne 
involves the following c_.i.er 


choice of the rotation ra ; 


* 00 — 003 Vc IT i 2D J — iS 


6 = 


x 


a . 

mm _ 


a 


2ca“ I A - 

Q it - ~ 


(35) 


* 3 r a >2 = 

e = 2 " l n J ^ 

2 5 max 


( 36 } 


• i •• j i_ ^ p'"ippc i "'", ess auanL — — ^ 3 ^ 

Here we have mtrocuced me c.m 

t drived bv the relations 

indeoendent of tne rota-ion . - — / a — a '"-~ 


<j, 3 = 2 a 5 2 


2 ( n . / p. ) 

v mm 7 max 


2 e/a max ) 


2a 6 g 2 Tb 2 (2 Ap /Si)' 


(37) 


The definition using $ t an( 3 


the factor 2, are for la^er 


convenience . 


The limit n min correspo 


r.ds to the strict requirement 


5 < 1, 


(38) 


so 


that the spiral center is inside the rotating react 


or. The 


limit 0 correspnds to the loose requirement 

max 
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z < 1 , 


„ _ 1 


soirallm 


so that the cenuriruga 
oeriod T. The radios changes 


is not excessive curing the 

£ • ■ _ i_ V* ! -» O 

c a factor e curing ^ 


- => v 3 i 


v US ■ 




_ ^ : - ; 0 Lz t the c^rt i-i-s 

oarticle density excess Ac- ^ ~ ' 

, _ r -accicn F cf the react:: 

spiral cucwarc, an: 

. ___ ^ wh ' ~h carricies starting there will ncc 

cross section a_ = a --- Wi1 a 

• . t p or necative lr, 

1 ■ . .3T v G- - • * JL - 

h ’’ t tne reactor wail - a ^ 

* *- a v ' t ^ ^ tne i r s c >- icn F •— - 

spiral inwarcs, anc *e 

■ — = - — '"nc oartiCj.es w* . . x — n ha.— 


the particles 
the reactor cross s 
not hit th 


:n a: 


. . — a i T • 1 IS mU 1 --- 

wall in tne prec=--n 3 

■ ^ ^<en as 10 5 sec (or 27.3 

more than a rotation pe-xOw a..- -- 

. vi _ _ - rn a - ‘his analysis makes nc use 

hours) for our examp_es. Nl.. - na - 

’ 3 o s e in:eoca ! Ley (3 3 ' . 


■ ^ T(— -j i—> r*n o w ~ p n its esn-.r: 

In both cases, tne area mvox/-- - s a 

the soiral center to tn- 


it the soiral center. The distance rrom w. 


- . -j n n • /K- I V I) T'r^o ^a 0 1 U 5 O E tflS Cl 

cylinder wall is ( - i x o ! * ~* iw ~ 


ircle is there: 


r F = (b - x o J )« 


- £ 


(40) 


The fraction r i = ^ ne A ^ - 
sma 1 last circle in Figure 3) t, 
largest circle). Thus 


o : 

the area 

of t h 

the 

cy l i nder 

cross 


section area Ct.ne 


F = 


irr F v -2e 

— L_ = (1 - 6) e 


77b' 


3 , - 2 . 


Z / J / - *• \ 

= ( 1 - 6 ) exp (- ♦ / & * 


(41) 
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15 


m tms 


1 3 S t 0XOreSSlOn -O j 


regains cons 


:ar.t while P. and 5 


vary . F is 

maximised 

by differentiation witn 

setting the 

derivative to zerc. Her.e = 

.3 



0 

II 

-o 


1 - a 


This is a c 

ubic equation in 6, warn a s_ngj.e 

allowed int 

erval 0 < 

5 < 1, for all positive 

Approximate 

analytic 

so lu t ions for <5 can be c 

either very 

small or 

very large; otherwise nu 

methods can 

be used, 

or o can be obtained usi 


r i ■ 
\ ^ ■ 


7 *5 1 1 


or crao'nical 


described below. 

Once the optimum 6 and F have been found, --- 2 p a .ui-J-a. 
set of rotating reactor parameters determining tr.e o values 
(37), the optimum rotation rate is given by equation (35), in the 

form 


a = n . / o 

mm 7 


(43) 


Using equation (42), wirn tne ce-in- w _ons 
be written as the alternative form 


(35) and (37), this cs 


.n 


si = a (l - 5) 1/3 / 

S 


(44) 


where 


rvM'h-' 1 r: h.Gl ? 

OF FOO?? QUA'JT 
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tr > 



°s - “mi:/* * hb ? 1 


( 45 ) 


/.. 0 ^ r . ic1 3 s ; :e and density, and of the fluid 
is independent or tne uar>-ic-_ 

_ _ -[ ^ W .-V h c = 10 3 cm/ sec 2 , b = 3 cm, and 

OUT 0 K bilup 1 3 r * - 1 - 

-i v- i.' ! 22 rpm. The expression (4s) 

-1/3 - 

is useful, because (1-5) ' c: 


orcoerties. For 


T = 10 5 sec, H s is 0.113 sec * oi 

^-1/3 


:s not increase s--*.-- 


from unity until 5 approac.nes ----i* 

e c-s’i *, n . is much less t.nar. 

For the easy case of ve.v s*.u r in 


. The aDDroximate an 

max 


alvtic solution is 


6 = <t> 


£ - <}>/2 , 


F = 1 — 3 4> 


(46) 


The correspon'ding optimum 


rotation rate 


for small $ is 


a = n ( 1 + <j>/3 ) , 


(47) 


or approximately • 1 ^ ie r ° u 


ation rate is independent of the 


>article and fluid prope 


rties so long as they ensure that 


is 


small • 

• .u J _ j 4 n ,.k e o-evious section leads to a 
The example introduced m 

very small * value. Using those values. 


ornoimi page is 

OF POO* DUALITY 
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n 


-4 "I 

= 7.4 x 10 sec 


mm 


^max 


-1 

= 2.U sec 


, 3 = 2 < a ,r,ih /! W 


2 ^ 4 x 10 


-7 


5 = <f> 


= 0.00625 , 


= ♦/: 


= 0.00312 , 


= 1 — 3 <t> 


= 0.9813 , 


rotation 


The oot iraum 

and is almost exactly - s 


rate is s' 


"jrately given by equc 


fl . The so 


iral center is ar 


I x I = bS = 0.019 
l o 


CTtl 


., bl » d<sol»cenent from the react 
, ich is an impercept.bl- ; - - f3r 


t-’no ni-her hand, the 




6 = 1 - 4 ’ • 

e = <4 3 /2 r 
p = exp (“4> 3 ) * 


( 43) 

% 

t ion ( 47 ) , 
(49) 

or axis . 
large <}> 

( 50 ) 
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„ , w i th a • much great 

This is the hard case or large — s, -- mm 


a r 


than ft 


'max 


use I ess. 


hard case of 

large part: 

The optimum a 

is only s. 

>s centrifuge 

so raoidly 

.y small, and 

the Rotary 

1 0 r m ediate vs — 

- , 3 

ues or <? , 

"able 1 may be 

pref erabl 


'min 


v Reactor is in 


O ~ ' 


to solving the 


, . - ; n egpef t ive to <5 near the 

since F -3 mss—* 

equation (42), Dc- 

, u ; i ttp^ups of 5 soaced 

«, tu p v a h'° shows optimal va-ue^ ui. 

optimum vaiae. xne ta_- 

, i-he corresocndmg values ol 

unif ormly from zero to one, w t. 


*3 = 6 3 / ( 1 " 6 ) 


e = $^/2 5“ , 


.2 -2 c 

F = ( 1 ~ <5 ) e r 


(1 - 5 ) 


- 1/3 


, , u > h t = nd interoolat ion , to determine 

The table can be used, wi m h-na me . 

• „ _ - 1_ ; - -ran notation 

values for. the other parameters, assuming a . * . 

rate, for any value of <t> (° r °f e or r ^‘ 

To illustrate the application of Table 1, consider the 

example given in the previous section for which the radius is 
calculated using = « mi n* ThiS fliV6S SUCCSSSi ° n 
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_ T 


a = a 

rricL 


= 1.42x10 cm < 


= n = 2 1/3 n s = 0.1454 sac 
max mm 5 


-i 


, 3 _ ^ 

9 - - / 


5 = 0 • / / 


e = ^ 3 /2 o 2 = 1.63 , 


F = (1 ~ 5) 


2 exz (_2e) = 0 . 0C 13 


a = n . /s = a_/U 

mi IV s 


- 0 


1/3 _ 3,19 39 se* 


-i 


• . « ,3 , int°’'colat ion in the table, and 

Here 5 is obtained ^rcxi <j» by m - . 


and F are then obtained from tne 


he fraction F of the 


. . a-ceotably snail, even at the 

particles not hitting tne v,ai- 

optimum rotation rac_us. 

.• f nr w'ni’ch a fraction F of 
To determine the largest radius for w... 

i the table, and obtain 

0.5 can be obtained, we again loo - 

approximately. 


5 — 0.2 , 


( j , 3 = 0.01 
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From equation (3/), 


(a/a„ J 6 - * 3 /2 - (0.4135) 


which g 


ives a radius limit ci 


diameter ci 


>tation rate is 


q /(I - 6) 1,/3 = 0.1278 sec 
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5. CONCLUSIONS 


, . , _ r-^t = tina about its 

i* <- a r'7 1 ■ naer oc iiu-j = 

A rotary reaczor is a ^-- iIUC - 

t)2 r'-->ies in suspension. Sucn a 

horizontal ax_s — n or- — 

i,,„ 3r . ir.i orosoheres in tne 

reactor has been used to procure 1 - 

. r- w -‘ -h results approaching those 

laboratory without: f icccula _--r. , w._n 

obtained in s p a c s . 

e „-,v n fallowed bv a particle in a 

4n ao;u:1 . 3 solution for tne pa^n 

, . j 'T’hp ns-h deoencs on whether 

rotating reactor has seen ooia.ne ■. 

the particle is heavier or lighter than the fluid. For 

counterclockwise rotation, a heavy particle spirals outwards 

around a spiral center displaced to the right or the axis. A. 

Ucht particle spirals inwards around a center displaced to the 

left. The spiral rotation rate is the same as the rotation rate 

of the reactor. The relative change in the ractus, for ea-h 

spiral rotation, is very small. 

Physically, the horizontal displacement of the spiral center 

from the axis of rotation is determined by the condition tnat tne 

, i - i i -, n( ~ 0 the viscous drag from the 

net weight or buoyancy shoulc w,alanc~ tne v- 

. . qnirallina outwards or inwards 

fluid flow past the particle, .xhe spiralling o 

_ / .■ q 4-Kp c^ntiri-fuccLl foircs cn 
is due to the centrifugal buoyancy (-• 

j • *.•«. *-v. e c°ntrifuaal force on the fluid 
the particle as compared wU tne c~. 

it displaces ) . 

There are two constraints on the rotation rate of the 
reactor. It must exceed the minimum value (26), to keep the 
spiral center inside the reactor. And it MS. not be muc g 

than the maximum value (33), since either too many particles hit 
the reactor wall or the particle concentration becomes excessive 
(for light particles spiralling inwards). 


r/V; 




u i rrv 


iw 


For both cases , 


there is a natural optimization problem, to 

^ ^ ^ the fraction c£ the reactor 

choose the rotation rate to 

. . . . , wh* : ch either will not 

cross section area wmcn ccn---.^ - 

- - * , , — ^ v -a pvnpr ' ^s 1 '- time , cir 

soiral cut to the reactor wa_- = ' 

have spiralled in without hitting me wa_i. 

*-n u, 0 c^t’Tiizetion Drool 0... +• 2, s - - - • - ^ ^ 

rorrestonc ing fraction v 


^ 0 1 e im 1 pi e l10 


optimum rotation rate and tne 
any set of parameters. Remarkably, the optimum rotation rate is 

independent of the particle and fit it properties, mr Sn^-1 
particles, and increases only a litt.e as the P c -'- iC - e 
increases . 

Table 2 presents the parameter values for the tnr~~ exa*,^ — s 


used in the texl 


The fluid properties, density c e r - n c , 


ate in each 


reactor diameter, and experimer.ua! a the 

case. Particle diameters in micrometers of 20, 234, and 117 are 

considered. The computed rotation rates are presented in rpm 
units for engineering convenience. The optimum rotation rate 
maximized the success fraction F of tne reactor cross section 
area. The spiral center displacements (from the reactor axis) 
are shown. The distance of each particle from the spiral c_nt_- 
changes by the indicated factor during the experimental time. 

For the first column, the particle diameter was chosen as a 
typical value of interest. The minimum rotation ra^e is much 
less than the nominal maximum, and a high success fraction is 
obtained, with a small spiral center displacement and a radius 
change factor close to unity. 
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, • ^ ^ v» 3 S cr.csen to 

In the second column, tne p^r 

• • m a »- i ~n *-at°s eaual. The optimum 

make the maximum anc minimum — 


is near tne 

. _ . r-r nf five, and 

. , ^ „ . , - V-> p >-, — o c; hv — £ a. w t. -«■ L o-L w * - 

] ] f-nc. n _ - uha.-o — — *-'j: 


rotation rate is somewhat higher, the spiral can 


rssctor wsll; tne ~ 


the sue 


- » rv 


' - verv 


. - /- 


■he oartiele c: 


or a 


The other ? 


chosen to 
er values 


For the third c: 
give a success rrac^ 
are shown . 

T ■ *. „ r - c * : ~ * o u - r a 1 v acoropr are - - 1 ■■» - - x 

The t ime T m tr.tr **0 i e -- * * 

microsphere processing. Note that the optimum - r — 1 

equation (44, and (45 varies only with the one-miri power of g/ = T 

and is practically independent of the particle rmiu= and 

provided (37) is small. Thus the optimum rotation rate is close t. 

rpm for latex microspheres, and is probably bet-e.., --- r P 

. . . , „ . « n .. r o ~ m : ■-* r c ? r a v i t y s i m u i a c - c . r . 

r n m ro* a very * - - * — 1 * r ■— ■** ■ 

Further work is required for a full application of this 
analysis to the production of mcnocisperse latex mi -- osphe.es . 

As described in Section 1, the monomer-swollen seed particles are 
buoyant in the early stages, anc converge on t..e _^__-l cen< " 

In the later stages, the particles shrink slightly anc become 
heavier than water, and the spiral center crosses the axis. The 
heavier particles centrifuge out as the polymerization approaches 

completion . 

The risk 'and extent of coagulum formation curing this 
process, and the deviations from a monodisperse s_^- v.isuribut’ 
due to variations in. the particle distribution, a.e unknown, 
can be expected that even if the particles collice near the 
spiral center as they centrifuge in during the early stages, the 


forces between the particles will be much smaller than for ^Lll 
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This is 


particles creamed to the top m a non rotating rea^-or. 
because the centrifugal accelerations are ex-remel/ sma„. Thus 


coaaulum 


formation might be minimal. In the late stages, as 


winn 


particles centrifuge out and hit the wall, there will be rurther 
particle collisions, with again the possibility of coagulum 
formation or sice dispersion. 

A further issue involves the beginning of procucuion, and 
the dynamics and thermodynamics of raising tne temperature, 
or without stirring. The stirring could £-SO be continued in<_o 
the polymerization stage. 

Finally, we assumed that tne so. Id-body rotation or tne fluid I £ 
not disturbed by the presence of the particles. However, recent 
experiments by Korn f eld show significant secondary flows driven by 
concentration variations, for C.3% suspensions of 50 micron latex 
particles, at rotation rates below 1 rpm. The dependence of this 
phenomenon on concentration is weak, while its dependence on rotation 
rate is very abrupt; at 1.4 rpm there is no observable modification 
from uniform rotation, while at 0.7 rpm the flow field and particle 
distribution are totally different. We plan further study of this 
effect. 
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Table 1. Parameter Values Corresponding to Uniformly S paced 6 Value s 
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Table 2 . Typical Parametric Values bo Illustrate Rotation Rale Liiail-? 



figure captions 


1 . A orototype Rotary Reactory to minimise c: 


ravi tat ional 


settling 


(or crear.inc) for the crowing -atax 


part icles . 


2. Schematic view a_or 

Rotary Reactor showing 


Ho - V ’ 


of symmetry c: me 
he coordinate sysren used 


par i;cias a “ A are stationary, supported a g a i n s t g r a / i : 
by the counterclockwise or clockwise flow. Other parti : i as 
move around A in circles, spiralling slowly outward or 
inward. Particles in the crescent hit the wall on the 
first cycle. Heavy particles initially between the 
concentric circles spiral out to hit the wall in time 7 . 
Heavy particles in the inner circle spiral out to fill the 
outer circle. Light particles initially filling the outer 
circle fill the inner circle after time T. In both cases 
we choose the rotation rate to maximize the area of the 
inner circle. 
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